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ABSTRACT 


This  report  presents  a  combined  analytical  and  experimental  investigation  of 
turbulent  heat  transfer  on  basic  and  composite  configurations  at  hypersonic  speeds. 
The  analytical  results  are  presented  in  Volume  I,  the  experimental  results,  including 
data- theory  comparisons,  are  presented  in  Volume  II,  and  computer  programs 
incorporating  the  analytical  methods  described  herein  are  presented  in  Volume  III. 

Two  analytical  approaches  are  presented:  the  prpr  method  and  the  turbulent 
nonsimilar  boundary  layer  method. 

The  prpr  method,  which  is  derived  from  the  boundary-layer  momentum  and 
energy  integral  equations,  is  recommended  for  predicting  turbulent  hea  ing  rates. 
Effects  of  dissociation,  pressure  and  wall  temperature  gradients,  three  -Himonsional 
flow,  and  nose  bluntness  are  included.  Simplified  methods  for  making  turbulent 
heating  estimates  using  a  slide  rule  or  desk  calculator  are  also  presented.  The 
computation  of  heating  rates  on  a  typical  reentry  configuration  in  flight  and  the 
extrapolation  of  test  data  from  ground  test  facilities  to  flight  are  described. 

The  turbulent  nonsimilar  boundary  layer  approach  offers  several  advantages 
over  previous  methods,  and  is  recommended  for  specific  parametric  studies  of 
turbulent  flows.  Calculations  made  using  this  method  have  been  restricted  to  ideal 
gases.  Modifications  can  be  made  to  include  real -gas  effects. 
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SECTION  I 


INTRODUCTION 


1.  REQUIREMENT  FOR  HEAT  TRANSFER  PREDICTION  METHODS 

The  need  for  better  methods  for  estimating  aerodynamic  heating  rates  is  contin¬ 
ually  increasing.  Relatively  crude  approximations  have  usually  been  adequate  for 
estimating  surface  temperatures  on  supersonic  aircraft  and  ablation  protection  re¬ 
quirements  for  ballistic  reentry  capsules.  However,  the  design  and  operating  limits 
of  hypersonic  cruise  and  maneuverable  reentry  vehicles  of  the  future  will  be  more 
sensitive  to  aerodynamic  heating  considerations.  This  will  require  more  accurate 
heat  transfer  prediction  methods. 

The  prediction  of  turbulent  heaiing  presents  a  particularly  serious  problem  for 
two  reasons.  First,  turbulent  heating  rates  are  usually  much  higher  than  the  corre¬ 
sponding  laminar  values.  Since  the  conditions  required  for  transition  from  laminar  to 
turbulent  flow  cannot  yet  be  predicted  with  confidence,  design  estimates  in  regions  of 
uncertainty  must  be  based  on  the  higher  turbulent  estimates.  Secondly,  the  complexity 
of  turbulent  flows  necessitates  the  use  of  simplifying  approximations  and  empiricism 
iu  formulating  the  fundamental  flow  equations;  thus,  methods  for  predicting  turbulent 
heating  rates  are  inherently  less  reliable  than  laminar  methods. 


2.  PRESENT  INVESTIGATION 

A  two-year  investigation,  including  both  analytical  and  experimental  studies,  was 
conducted  to  provide  and  verify  methods  for  predicting  turbulent  heating  rates  on 
basic  shapes  and  composite  bodies  at  hypersonic  speeds.  Results  of  the  analytical 
studies,  including  recommended  methods,  are  presented  in  Volume  1  of  this  report. 
Experimental  results  and  data-theory  comparisons  are  presented  in  Volume  II. 

Two  basic  methods  are  presented.  The  first,  the  PrPr  program,  is  recommended 
for  making  turbulent  heating  estimates.  The  second,  the  turbulent  nonsimilar  program, 
represents  a  new  approacn  in  treating  turbulent  flows,  and  is  intended  for  basic  studies 
of  turbulent  boundary  layer  phenomena.  Both  have  been  programmed  for  the  SRU  1108 
and  IBM  7094  digital  computers.  A  description  of  these  computer  programs  is  given 
in  Volume  III  of  this  report. 

a.  PrPr  Method 


The  PrPr  method  for  predicting  heat  transfer  and  skin  friction  for  both  laminar 
and  turbulent  flows  was  developed  by  Richard  A.  Hanks  of  The  Boeing  Company  in  the 
course  of  the  X-20  program.  Modifications  and  refinements  of  this  method  were  made 
under  a  subsequent  NASA  contract  (Reference  1).  The  method  presented  here  is 


1 


essentially  the  same  as  that  presented  in  Reference  1,  although  some  modifications 
to  the  computer  program  were  made  during  the  present  study. 

The  derivations  of  the  basic  p  .pr  equations  presented  in  References  1  and  2 
were  based  on  solutions  of  the  boundary- layer  momentum  integral  equation.  This 
equation  was  transformed  into  an  equivalent  incompressible  form  using  a  coordinate 
transformation  suggested  by  Mager  in  Reference  3.  The  transformed  equation  was 
solved  to  obtain  an  expression  for  skin  friction.  The  corresponding  heat  transfer 
equation  was  then  obtained  using  a  generalized  form  of  the  Reynolds  analogy. 

A  new  derivation  leading  to  the  same  result  for  heat  transfer,  but  based  on  a 
solution  to  the  transformed  boundary -layer  energy  integral  equation,  is  presented  in 
Appendix  A.  The  heat  transfer  equation  obtained  from  the  solution  of  the  energy 
integral  contains  boundary-layer  thickness  parameters  and  reference  densicy  and 
viscosity  terms  as  undefined  functions.  These  functions  represent  effects  of  fluid 
property  variations,  finite  stream  wise  and  crossflow  pressure  gradients,  and  stream- 
iine  divergence.  For  laminar  flow,  these  function  vere  evaluated  by  exact  solutions 
for  self-similar  boundary  layers.  Analytic  expressions  were  found  that  agree  with 
essentially  all  of  the  exact  similar  solutions. 

The  derivation  of  the  basic  prpr  heat  transfer  equation  and  the  correlations  used 
in  defining  the  undefined  functions  for  laminar  flow  are  presented  in  Appendix  A. 
Expressions  for  evaluating  the  turbulent  functions  are  presented  in  the  following 
sections. 

Because  of  the  overall  complexity  of  the  PrPr  equations,  this  method  is  not 
recommended  for  making  hand  calculations.  Handbook  methods  for  estimating  both 
laminar  and  turbulent  heating  rates  using  the  prpr  method  are  presented  in 
References  4  and  5.  1  Simplified  approximations  to  the  prpr  equations  permitting 
hand  calculations  arc  presented  in  Appendix  B  of  this  report.  Predictions  obtained 
using  the  simplified  equations  are  usually  within  5%  of  the  computer  results. 

b.  Turbulent  Nonsinular  Method 


The  turbulent  nonsimilar  method  is  an  extension  of  the  laminar  nonsimilar  method 
reported  in  References  1  and  (>.  In  this  method  a  semi-empirical  expression  is  used 
to  establish  the  Reynolds  stress  term  appearing  in  the  momentum  equation.  The 
corresponding  conduction  term  in  the  energy  equation  is  then  related  to  the  Reynolds 
stress  by  an  effective  turbulent  Prandtl  number. 

The  turbulent  nonsimilar  method  for  flows  was  developed  before  the  present  study 
was  started.  During  this  investigation  it  was  extended  to  include  compressibility  effects. 


1  Reference  4  supplements  Reference  5.  Except  for  delta  wings  at  angle  of  attack, 
turbulent  heating  estimates  obtained  from  these  reports  are  nearly  identical. 


The  analysis  was  conducted  to  provide  a  method  for  computing  velocity  and  enthalpy 
profiles  in  a  compressible  turbulent  boundary  layer  including  the  effects  of  pressure 
and  wall  temperature  gradients. 

The  derivation  of  the  turbulent  nonsimilar  equations  is  presented  in  Appendix  C  of 
this  report. 


SECTION  II 


BASIC  METHODS 


In  contrast  to  the  well-developed  theoretical  methods  available  for  laminar  flows, 
no  truly  analytic  basis  exists  for  treating  turbulent  boundary  layers.  The  unsteady 
velocity  fluctuations  greatly  complicate  the  mathematical  description;  consequently, 
ali  turbulent  methods  are  founded,  at  least  in  part,  on  empirical  correlations. 


1.  GENERAL  APPROACHES 


Most  turbulent  analyses  are  based  on  the  assumption  that  turbulent  flows  are 
analogous  to  laminar  flows  with  special  viscous  and  heat  conduction  properties.  These 
special  properties  can  be  expressed  in  either  the  fundamental  equations  of  motion  in 
differential  form  or  in  the  boundary  layer  momentum  and  energy  integral  equations. 

An  example  of  the  first  type  is  furnished  by  the  well-known  mixing  length  theory. 

a.  Mixing  Length  Theory 


For  incompressible  flow,  Prandtl  hypothesized  that: 


rT  =  pr 
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The  mixing  length  is  denoted  by  I .  Prandtl  further  suggested  that  l  was  nearly 
proportional  to  the  distance  from  the  wall.  Equation  (1),  using  the  approximation  that 
S.  is  proportion  to  y,  was  found  to  provide  good  incompressible  skin-friction  estimates 
and  velocity  profiles  when  the  shear  was  assumed  to  be  constant  with  y  (Reference  7). 
This  latter  assumption  is  clearly  unrealistic  for  exterior  boundary  layers,  since  the 
Reynolds  stress  must  approach  zero  at  the  edge  of  the  boundary  layer. 


Many  other  methods  using  mixing  length  theory  have  been  reported,  including  a 
widely  used  method  for  compressible  flows  developed  by  van  Driest  (Reference  8). 
However,  most  of  these  methods  are  based  on  assumed  velocity  or  shear  profiles.  In 
addition,  restrictive  assumptions  regarding  flow  similarity  have  been  required  in  order 
to  solve  the  flow  equations,  thus  limiting  these  analyses  to  flows  with  no  pressure  or 
wall  temperature  gradients.  Fewer  assumptions  are  required  by  the  turbulent  non- 
similar  method,  which  is  described  in  Appendix  C.  The  method  utilizes  a  shear 
correlation  similar  to  the  mixing  length  correlation  given  by  Equation  (1),  but  removes 
many  of  the  restrictions  imposed  by  previous  analyses.  By  solving  the  boundary- layer 
conservation  equations  in  partial  differential  form,  velocity  and  temperature  profiles 
are  then  obtained  as  part  of  the  solution.  No  limitations  arc  imposed  regarding  pres¬ 
sure  or  wall  temperature  gradients  except  that  separated  flows  cannot  be  treated. 
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b.  Integral  Methods 


Boundary  layer  methods  derived  from  momentum  and  energy  integral  equations 
are  often  much  easier  to  use  than  methods  requiring  the  exact  solution  of  the  flow 
equations.  The  energy  integral  equation,  Eq.  (A-9),  derived  in  Appendix  A  of  this 
report  and  the  momentum  integral  equation  derived  in  Appendix  B  of  Reference  2  can 
be  expressed  in  the  forms  given  below: 
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The  equations  given  above  are  exact,  except  for  the  usual  boundary  layer  assumptions, 
and  are  valid  for  both  laminar  and  turbulent  flows.  However,  the  solutions  to 
Equation  (2)  are  quite  different  for  laminar  and  turbulent  flows,  since  the  velocity  and 
enthalpy  profile  parameters  used  in  solving  the  equations  are  different. 


2.  PrPr  METHOD 
a.  Basic  Formulation 


The  PrPr  equations  are  derived  in  Appendix  A  from  the  boundary  layer  energy 
and  momentum  integral  equations.  The  basic  correlations  are  of  Blasius  type,  and 
are  given  by: 
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The  symbols  in.  Cm,  and  Cq  are  profile  parameters,  6  is  the  momentum  thickness, 

and  Q  is  the  energy  thickness.  Momentum  and  energy  thicknesses  are  defined  in 

Appendix  A,  as  are  methods  for  evaluating  the  reference  density-viscosity  product 

P*  P 
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The  integral  equations  are  reduced  to  an  incompressible  form  using  a  modified 
Stewartson  transformation  and  are  given  in  detail  in  Appendix  A.  The  resulting 
equations,  in  physical  coordinates,  for  II,  rw/ue,  and  6  are: 
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where  the  equ. valent  distance  parameters  xecj,  and  Secj  are  defined  by: 
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The  equivalent  distance  xec,  T  is  used  only  in  computing  heat  transfer,  and  S^x 
only  in  computing  skin-friction  and  momentum  thickness. 


Methods  for  evaluating  the  terms  appe  ring  in  Equations  (5)  through  (9)  for 
laminar  flow  are  given  in  Appendix  A.  Expressions  for  evaluating  the  equations 
given  in  the  following  sections  reflect  the  latest  state-of-the-art,  and  are  essentially 
the  sam  :  as  reported  in  References  1  and  4. 

For  turbulent  flow,  the  profile  constants  Cx  and  m  are  discussed  later  in  this 
section.  Methods  for  evaluating  the  other  parameters  are  presented  in  the  following 


sections. 

Symbol  Definition  Section 

p  n  reference  density-viscosity  product  III 

p  reference  stagnation  viscosity  III 

influence  of  Prandtl  number  on  heat  transfer,  q/q  ,  III 

Pr  a  =  1 

£  influence  of  atomic  diffusion  on  heat  transfer,  III 

q/qN 

‘  Le  =  T 

J  profile  parameter  reflecting  pressure  gradient  effects  IV 

on  heat  transfer 

P  profile  parameter  reflecting  pressure  gradient  effects  IV 

on  wall  shear  stress 

r  streamline  divergence  parameter  due  to  body  geometry  V 

f  streamline  divergence  parameter  due  to  crossflow  V 

pressure  gradients 

E  crossflow  momentum  thickness  V 

b.  Evaluation  of  Cx  and  m 


The  profile  parameters  Cx  and  m  are  by  definition;  as  in  the  laminar  case, 
independent  of  pressure  gradients  and  flow  compressibility;  hence,  these  parameters 
can  be  evaluated  on  the  basis  of  incompressible  flat-plate  flow.  For  this  case 
pr  pr  =  pe  pe,  pQ  -  pg!  anc'  xecj  =  S  =  x;  thus,  Equations  (5)  and  (6)  simplify  to: 
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When  Cx  =  .0290  and  m  =  4,  Equation  (10)  corresponds  to  the  Colburn  equation,  and 
Equation  (11)  becomes  the  familiar  Blasius  expression.  No  heat  transfer  data  are 
available  for  comparative  purposes,  since  the  very  low  heating  rates  associated  with 
nearly  incompressible  gas  flows  cannot  be  accurately  measured.  Comparisons  with 
skin  friction  data  show  that  Cx  and  m  are  themselves  functions  of  the  Reynolds 
number.  To  include  this  variation  with  Reynolds  numbers  would  greatly  complicate 
computations;  hence,  an  alternate  formulation  for  the  skin  friction  expression  was 
sought. 

After  a  survey  of  several  proposed  incompressible  friction  formulas,  a  minor 
modification  of  the  Schultz-Grunow  equation  (Reference  9;  was  selected: 
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The  modification  is  the  addition  of  the  constant  3000  to  the  Reynolds  number.  This 
modification  was  made  in  order  to  provide  more  realistic  values  of  Cf  at  Reynolds 
numbers  below  10")  (see  Figure  1). 

The  modified  and  unmodified  expressions  are  shown  in  Figure  1,  together  with 
other  available  methods.  As  shown,  there  is  little  difference  between  the  various 
methods,  except  that  the  Blasius  equation  falls  low  at  high  values  of  Reynolds  number. 
Equation  (12)  was  originally  selected  because  of  its  slight  conservatism,  although  any 
of  the  other  expressions  could  have  been  used. 

The  form  of  Equation  (12)  does  not  lend  itself  to  calculations  in  the  framework  of 
Equations  (10)  and  (11)  due  to  the  variation  of  m  with  Reynolds  number.  However, 
comparisons  have  been  made  that  show  that  nr  =  4  is  an  adequate  approximation  for 
eraluating  geometric  effects.  For  example,  if  m  is  evaluated  at  particular  values 
of  Reynolds  number  using  Equation  VT2)  the  following  comparisons  are  obtained: 
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Supporting  experimental  data 
are  presented  in  Reference  7 


COMPARISON  OF  INCOMPRESSIBLE  TURBULENT  SKIN-FRICTION  FORMULAS 


Thus  the  effect  of  variations  in  ni  is  seen  to  be  small.  Accordingly,  m  =  4.  0  was 
selected  for  calculation  of  geometric  effects  presented  in  Sections  III,  IV,  and  V. 

However,  since  sizeable  errors  can  result  from  using  m  =  4  in  Equations  (5)  .\nd 
(6),  the  basic  form  of  Equation  (12)  has  been  adopted.  The  new  expressions  replacing 
Equations  (5)  and  (<>)  are: 
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The  definitions  of  xeC(  j  and  Se(1(x  are  given  by  Equations  (8)  and  (9),  respectively. 

The  Schultz-Grunow  expression  for  momentum  thickness  of  an  incompressible 
flat-plate  boundary  layer  is  given  in  Reference  9  by: 
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SECTION  m 


COMPRESSIBILITY  AND  DISSOC LA T ION  EFFECTS 


The  influence  of  compressibility  and  dissociation  effects  on  heat  transfer  and 
skin  friction  is  usually  accounted  for  by  basing  the  density  and  transport  properties 
(i.e.  ,  coefficients  of  conductivity  and  viscosity)  on  some  reference  thermal  param¬ 
eter,  For  example,  reference  temperatures  (or  enthalpies)  based  on  a  weighted 
average  of  the  wall,  adiabatic  wall,  and  ooundary  layer  edge  temperatures  (or 
enthalpies)  are  commonly  used  as  a  thermal  reference.  The  PrMr  method  differs 
from  the  common  reference  temperature  methods  in  two  respects.  First,  a  reference 
density -viscosity  product  PrPr  is  evaluated  instead  of  a  reference  temperature  and 
secondly,  a  reference  stagnation  viscosity  pQ  must  be  evaluated.  Other  correction 
terms  appearing  in  Equation  (13)  reflecting  effects  of  compressibility  and  dissocia¬ 
tion  are  expressed  by  the  Prandtl  number  function  Fpr  and  the  Lewis  number 
function  £. 


1.  EVALUATION  OF  REFERENCE  DENSITY -VISCOSITY  PRODUCT 

The  reference  density-viscosity  product  for  turbulent  flow  is  taken  to  be  the 
laminar  value  defined  by  Equation  (A -42)  in  Appendix  A  and  plotted  in  Figure  2  (for 
air  in  chemical  equilibrium,  the  density-viscosity  ratios  required  to  obtain  prMr 
can  be  obtained  from  Figure  3).  This  basic  identity  is  suggested  by  the  fact  that 
PrPr  appears  only  in  connection  with  the  laminar  shear  terms  of  the  turbulent 
boundary  layer  equations^.  Further  justification  is  provided  by  the  excellent  agree¬ 
ment  between  estimates  from  the  PrPr  method  and  most  experimental  data  obtained 
from  several  facilities  covering  a  wide  range  of  test  conditions. 

Comparisons  of  prpr  predictions  v\ith  flat  plate  heat  transfer  data  obtained  from 
the  Cornell  Aeronautical  Laboratory  (CAL)  48"  shock  tunnel  and  reported  in  Reference 
1  are  shown  in  Figure  4.  The  PrMr  correlations  given  by  Equation  (A -42)  provide 
good  agreement  with  X-15  flight  data  as  shown  in  Figure  5.  The  X-15  comparisons 
are  especially  noteworthy  because  of  the  lack  of  agreement  with  other  widely  used 
methods.  These  compar'sons  demonstrate  good  agreement  between  theoretical  pre¬ 
dictions  and  experimental  measurements.  More  recent  flat  plate  data  obtained  from 
the  CAL  shock  tunnels  and  reported  in  Reference  10  are  shown  in  Volume  II  of  this 
report  to  be  substantially  higher  than  the  prpr  estimates.  The  reason  for  the 
discrepancy  between  the  two  sets  of  data  is  not  fully  understood  at  this  time.  Some 
of  the  possible  causes  are  discussed  in  Volume  II. 


1 1 


2  E.g.  ,  Equation  (13)  in  Reference  3. 


Figure  2:  REFERENCE  DENS  I TY -VIS COS  I TY  PRODUCT  PARAMETERS 
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Figure  3:  VISCOSITY  FUNCTION  FOR  AIR  IN  CHEMICAL  EQUILIBRIUM 
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Figure  4:  HEAT  TRANSFER  ON  A  FLAT  PLATE  IN  A  SHOCK  TUNNEL 
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TURBULENT  HEAT  TRANSFER  THEORY  COMPAR  ISONS  WITH  X-15  FLIGHT  TEST  DATA 


2.  REFERENCE  STAGNATION  VISCOSITY 


Since  the  reference  stagnation  viscosity  does  not  appear  in  the  laminar  equations, 
no  information  regarding  its  evaluation  can  be  obtained  by  examining  the  laminar 
solutions.  The  reference  stagnation  viscosity  is  assumed  by  Mager  (Reference  3)  to 
be  the  viscosity  evaluated  at  stagnation  conditions.  For  real-gas  How  with  the  vis¬ 
cosity  dependent  on  the  pressure  it  seems  more  realistic  to  consider  the  local  flow 
composition  rather  than  the  composition  corresponding  to  stagnation  conditions. 
Accordingly,  p0  is  calculated  with  the  Sutherland  law  and  pr  using  the  value  of 
specific  heat  corresponding  to  PrPr.  The  result  is: 


+  200 


(19) 


When  prpr  and  local  pressure  are  known,  the  compressibMity-temperature  product 
ZT  can  be  obtained  for  air  in  chemical  equilibrium  from  Figure  3.  The  corresponding 
reference  enthalpy  ir  can  then  be  determined  using  Figure  G. 


3.  EVALUATION  OF  PRANDTL  NUMBER  AND  LEWIS  NUMBER  PARAMETERS 

As  with  the  PrPr  correlations,  the  effects  of  Prandtl  number  on  turbulent  flat 
plate  heat  transfer  are  assumed  to  be  identical  to  the  laminar  case,  the  correction 
term  is: 


(20) 


Where  a  is  the  partial  Prandtl  number  for  translation,  rotation,  and  vibration.  The 
subscript  r  denotes  that  a  is  evaluated  at  the  enthalpy  and  composition  correspond¬ 
ing  to  PrPr.  For  equilibrium  air,  ct  car.  be  obtained  from  Figure  7. 

Similarly,  no  analytic  basis  has  been  found  for  determining  the  influence  of 
atomic  diffusion  on  turbulent  neating  rates.  Furthermore,  the  combination  of  high 
pressures  and  temperatures  required  to  obtain  significant  dissociation  of  air  in 
turbulent  boundary  layers  cannot  be  obtained  in  most  present-day  ground  test 
facilities. 


l(i 


COMPRESSIBILITY-TEMPERATURE  PRODUCT 


FOR  AIR 


In  the  absence  of  a  rigorous  method,  the  influence  of  dissociation  is  assumed  to 
be  the  same  for  turbulent  and  laminar  flows.  The  laminar  expression  for  £  is  given 
by  Equation  (A -41)  as: 


£  =  i  +  (n_  -52 

'  Le 


(21) 


A  plot  of  £  as  a  function  of  enthalpy  and  pressure  for  equilibrium  air  and  a  Lewis 
number  of  1.4  is  presented  in  Figure  20  in  Appendix  A. 
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SKCTIOX  IV 


STREAMWISE  PRESSURE  GRADIENT  EFFECTS 


Except  lor  boundary  layer  profiles,  the  effects  of  pressure  variations  on  heating 
and  skin  friction  are  reflected  in  the  product  P,.pr ue  appearing  in  the  integrand  in 
Equations  (S)  and  (9).  Changes  in  boundary  layer  profiles  resulting  from  pressure 
gradients  are  accounted  for  in  Jy  and  Py  appearing  in  the  same  equations. 

At  the  time  that  these  correlations  were  made  (1964  and  1965),  no  reliable  methods 
were  available  for  estimating  turbulent  profile  alterations  due  to  streamwise  pressure 
gradients.  Consequently,  the  selection  of  expressions  for  evaluating  Jy  and  P.j.  are 
based  on  comparisons  with  experi mental  data.  The  general  form  of  these  equations 
is  based  on  the  corresponding  laminar  relations  presented  in  Appendix  A.  The 
expressions  selected  are: 

J  =  JT  =  |l  -  .  718  ( -  F  F_  -1)1  when  0  >0  (22) 

1JL1  P  .  S  L  y  S  i  S 

and. 

vjl  <2S> 

Also; 

PT=PL=JL  '21> 

Expressions  for  evaluating  F^  g  and  F^,  ^  are  given  in  Appendix  A  by 
Equations  (A- 50)  through  (A- 53). 
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SECTION’  V 


THREE-DIMENSIONAL  EFFECTS 


Although  three-dimensional  effects  on  turbulent  heating  are  significantly  less  than 
for  laminar  flow,  these  effects  can  still  be  substantial.  In  particular,  sizeable  in¬ 
creases  in  heating  due  to  crossflow  pressure  gradients  can  be  experienced  on  leading 
edges  and  along  the  stagnation  line  of  axisvmmetric  bodies  at  angles  of  attack. 


In  the  present  analysis,  three-dimensional  effects  on  heat  transfer  are  reflected 
in  the  equivalent  distance  defined  by  Equation  (8).  The  influence  of  streamline  diver¬ 
gence  due  only  to  body  geometry  is  determined  by  r,  and  that  due  to  pressure  gradi¬ 
ents  normal  to  the  direction  of  flow  by  f.  The  combined  effects  of  body  geometry  and 
crossflow  pressure  gradients  on  the  streamline  divergence  at  the  boundary  layer  edge 
is  defined  by  A,  where: 


A  =  rf 


Streamlines 


x 


The  reason  for  treating  the  two  causes  of  three-dimensional  flow  separately  is 
that  crossflow  pressure  gradients  distort  the  crossflow  velocity  profile,  whereas  body 
geometry  does  not.  The  evaluation  of  the  influence  of  this  distortion  on  heat  transfer 
is  easier  when  crossflow  effects  are  considered  separately  from  geometric  effects. 


1.  GEOMETRIC  EFFECTS-  r  FACTOR 

The  streamline  divergence  parameter  r  is  generallyused  in  connection  with  axisym- 
metric  flows,  and  is  then  defined  as  the  circular  radius  of  a  body.  In  this  analysis  r  is 
considered  to  be  the  distance  between  two  adjacent  streamlines  at  the  edge  of  the  bound¬ 
ary  layer  over  the  respective  surface,  and  in  the  absence  of  crossflow  pressure 
gradients.  For  most  applications  r  can  also  be  considered  to  be  proportional  to  the 
body  radius  of  curvature  in  the  plane  normal  to  the  streamline. 

An  exception  to  this  rule,  delta  wings  at  angle  of  attack,  is  discussed  later  in  this 
section. 
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The  simplest  example  of  three-dimensional  flow  occurs  on  an  unyawed  cone  in  a 
uniform  hvpersonic  flow.  Noting  that  pressure  is  constant,  and  considering  that  the 
wall  temperature  is  also  constant,  the  product  prprue  appearing  in  Equations  (8)  and 
(9)  is  constant.  Also,  from  Se  on  IV  it  is  seen  that  Jj  and  Pj  are  unity;  hence. 
Equations  (s)  and  (9)  reduce  to: 


(25) 


The  streamline  divergence  parameter  r  is  proportional  to  the  circular  radius, 
and  is  therefore  proportional  to  x.  Considering  that  II  is  approximately  proportional 
to  xeq(  t~“,  it  is  seen  that: 


For  other  types  of  unyawed  axisvmmetric  bodies,  the  pressure  is  not  usually  constant, 
and  the  variation  in  prMrue  must  be  considered. 


2.  CROSSFLOW  PRESSURE  GRADIENTS  -f  FACTOR 

Estimating  the  effects  of  crossflow  pressure  gradients  is  more  difficult  than 
effects  of  body  geometry  because  of  rotation  of  the  streamline  patterns  in  the  boundary 
layer.  For  example,  consider  the  differences  in  the  streamline  pattern  on  an  unyawed 
cone  and  a  swept  cylinder  stagnation  line  sketched  below: 

Streamlines 

St  ream- 
lines  at 
all  values 
of  y 

Conical  Flow  Swept  Cylinder  Flow 

The  crossflow  streamline  divergence  parameter  f  can  be  considered  as  proportional 
to  the  distance  between  two  adjacent  streamlines  at  the  edge  of  the  boundary  layer  in 
the  absence  of  geometric  effects.  Methods  of  estimating  f  are  presented  later  in  this 
section. 


The  influence  of  the  streamline  rotation  illustrated  above  is  reflected  in  the  cross- 
flow  momentum  thickness  ratio  K.  When  streamlines  are  parallel  throughout  the 
depth  of  the  boundary  layer  at  a  given  station,  then  u/ue  =  v/ve.  It  is  seen  from 
Equation  (A-25)  that  E  is  unity  for  this  case,  and  the  effect  of  f  becomes  identical 
to  that  of  r  in  the  expression  for  heat  transfer.  Equation  (5). 

a.  Evaluation  of  Bp 

The  behavior  of  E  in  turbulent  flow  can  be  described  only  qualitatively,  and  most 
published  analyses  neglect  its  effect.  However,  its  effect  is  usually  to  increase  heat¬ 
ing  rates  and  is  therefore  included  in  the  present  method.  As  in  the  streamwise 
pressure  gradient  case,  the  turbulent  values  are  based  on  modifications  of  the  corre¬ 
sponding  laminar  correlations.  However,  unlike  the  streamwise  parameter  Jp,  E^ 
is  strongly  influenced  by  Mach  number,  so  that  a  dual  modification  is  required. 

Considering  incompressible  flow,  it  is  seen  from  the  definition  of  Equation  (A-25) 
that  the  upper  limit  on  E  is  6*/0  unless  the  crossflow  velocity  component  v  within 
the  boundary'  layer  exceeds  the  external  value.  Laminar  solutions  (Reference  11)  show 
that  velocity  overshoots  (v/ve  >  1)  do  not  occur  for  cold-wall  zero-Mach-number  flow, 
hence  a  correction  factor  of  the  following  form  is  suggested: 

E  (6*/6) 

T~  -  <6*78)7  =  c°"s,ant  <27> 

L,  o  L 

However,  an  inconsistency  is  noted  in  that  Ej  Q  should  be  unity  when  £p  0  is  unity, 
since  this  condition  implies  that  no  profile  distortion  occurs  due  to  crossflow'  pressure 
gradients  (i.e. ,  v/ve  =  u/ue).  This  inconsistency  is  avoided  by  adopting  the  following 
expression: 


-  1 


-  1 


0.77 


(28) 


The  constant  .77  was  selected  primarily  on  the  basis  of  experimental  heating  rates  on 
swept  cylinders. 

In  Equation  (28),  El  0  is  just  Ep  evaluated  for  Mach  number  equal  to  zero.  For 
Mach  number  zero.  Equation  (A -58)  reduces  to: 
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Assuming  that  F 


=  1.0,  and  using  Hquations  (28)  and  (A-49), 
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where 


and 


K_  =  1  -  .  55  (  v/T 
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\exp  K 


c,  o ' 
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c,  o 


(ZT)  /(ZT) 

m.  o  e,SL 


(30) 


c.xp  K  =  0  when  N  <  .  05  and  .99  <  NT  <  1.01, 

2 

exp  K  =  -  .  194  exp  ( -  -  X(X-l)]  when  .  05  <  X  <  .  99, 

O 

exp  K  =.  194  exp  (-^(X-l)]  when  X>  1.01,  and  X  =  (x/rf)  (9rf/9x). 

u 


The  effect  of  Mach  number  on  Kj  was  determined  from  observation  of  empirical 
trends  in  swept  cylinder  stagnation-line  turbulent  heat-transfer  data,  as: 


E. 


ET,o  \  EL,  o  , 


(31) 


The  final  expression  for  Ej  now  becomes: 


E  =  1  +  .  55  (  >/l  +  F  -  l)  ( 22  ) 

T  ’  2,  c,  o  '  ’•  c,o/ 


exp  K 


,1+-718  ('/lH  Fr,c.o-1)K,„rPK. 


(32) 


where  the  sign  on  the  exponents  is  plus  when  9(rf)/9x  is  positive,  and  negative 
when  9(rf)/9x  is  negative.  Although  the  analysis  leading  to  Equation  (32)  is  based 
on  order-of-magnitude  approximations,  heat  transfer  estimates  obtained  using  this 
method  are  shown  in  Figure  8  to  agree  reasonably  well  with  test  data  obtained  at  the 
stagnation  line  of  a  swept  cylinder. 
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Shock  funnel 


Figure  8:  CROSSaOW  PRESSURE  GRADIENT  EFFECT  ON  CYLINDER 
STAGNATION-LINE  TURBULENT  HEAT  TRANSFER 


b.  Evaluation  of  f 


1)  General  Equation  for  Axisymmetric  Bodies 

The  streamline  divergence  parameter  f  is  obtained  from  a  solution  of  the  cross- 
flow  momentum  equation  for  inviscid  flow  given  by: 


e  ,  e  e 

P  u  P  v 

e  e  9x  e  e  9y  9y 


(33) 


The  coordinate  system  is  defined  in  the  sketches  shown  below: 


Differentiating  Equation  (33)  with  respect  to  y,  and  noting  that  at  the  stagnation 
line  vg=  9ue/9y=  9pe/9y=  0,  then: 

9  9  2 

/3v  \'  a  v  a  p 

Pe  (ir)  4  ”eUe  3x8y  =  '  <3 

For  most  configurations,  pressures  are  more  easily  expressed  as  a  function  of 
x  and  0.  Equation  (34)  is  now  transformed  to  (x,  0)  coordinates.  The  relations  for 
this  transformation  are: 


90  =  1  30 

9y  x  r 


(35) 


90  j  =  90  _  0  dr  90 

9x  I  9x  r  dx  90 

I  y  0  x 

In  the  transformed  system,  Equation  (34)  becomes: 


9v  dr/dx 
_ e _ 

90  r 


(36) 


(37) 
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Equation  (37)  cannot  be  easily  solved  unless  it  is  assumed  that  32ve/3x30  =  0. 
This  assumption  is  clearly  valid  for  yawed  cones,  and  numerical  results  indicate  that 
errors  introduced  by  neglecting  this  term  on  yawed  ogives  are  small.  Accordingly, 
Equation  (37)  is  now  simplified  to: 
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where 

V  =  V  /u 
e  e 


-  ^  ,  2 
P  =  P  /  p  u 
e  e  e 


The  solution  to  Equation  (38)  is  given  by  the  familiar  quadratic  equation, 
resulting  in: 


3  v 
30 
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dr  ±  \ 
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30 2  . 

(39) 


The  proper  sign  in  Equation  (39)  is  positive,  as  seen  in  the  special  case  for  the  stag¬ 
nation  line  of  a  swept  cylinder  (3r/3x  =  0)  for  which  the  positive  value  of  3v/30  shows 
diverging  streamlines.  Thus: 


It  is  seen  from  Equation  (39)  that  when  32P/30^  =  0,  then  3v/3 0  =  3r/3x. 

Thus,  3v/30  reflects  the  streamline  divergence  due  to  body  geometry  as  well  as 
crossflow  pressure  gradients.  This  result  appears  to  be  inconsistent  with  Equation 
(A- 10),  which  shows  that  v  includes  only  crossflow  effects.  The  reason  for  this 
apparent  contradiction  is  that  in  the  analysis  shown  above,  v  is  defined  as  the  velocity 
component  normal  to  lines  of  constant  y,  but  in  Appendix  A,  v  is  the  component 
normal  to  lines  of  constant  0.  With  the  terminology  used  in  this  section: 


13r+l3f  =  13v 
r  3x  f  3x  r  30 


or. 


where 
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sin  h  =  3r/ 3x 


(40) 


2)  Evaluation  of  3^p/9p2 

Near  the  stagnation  line  the  pressure  variation  with  p  is  defined  by: 

2 

^  sin  6 

P  VT,SI.  2 

sln  6sl 

where  5  is  the  angle  of  the  suifac  :  with  respect  to  the  free-stream  flow.  For  an 
ax i symmetric  body: 

sin  6  =  cos  a  sin  7}  +  sin  a  cos  ij  cos  p 
Using  Equations  (41)  and  (42),  it  is  easily  shown  that: 


2  p  sin  a  cos  t) 
sin  (a  +  t)  ) 


Methods  for  evaluating  P  for  a  few  simple  geometries  are  presented  in  Section  VIII. 


3.  SPECIAL  CASES 

a.  Swept  Cylinder  Stagnation  Line 

For  the  stagnation  line  of  swept  cylinders  dr/dx  =  0,  thus  only  crossflow  effects 
need  be  considered.  For  this  case,  Equation  (40)  reduces  to: 

1  3f  1  31 2 3P 

f  3x  r  302  *44 

2-  2 

3  P/3p  is  found  from  Equation  (43).  Then, 


2  — 
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f  =  exp 


(45) 


For  this  case  PrPr ue  can  normally  be  considered  constant,  and  is  seen  to  be 
unity  { 0S  =  0).  Equation  (8)  now  simplifies  to: 
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For  very  large  values  of  x  (e.  g. ,  infinite  cylinder),  Equation  (46)  becomes: 
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b.  Yawed  Cone  Stagnation  Line 

The  flow  along  the  stagnatio .■>  line  of  a  yawed  cone  differs  from  that  on  the  swept 
cylinder  only  in  that: 


dr 

dx 


sin  rj 


where  rj  is  the  cone  half-angle  (degrees).  For  this  case,  Equation  (40)  becomes: 
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A  solution  for  f  is  given  by: 
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Again,  prjiru  can  be  considered  constant,  and  Jy  assumed  to  be  unitv.  Hence, 
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c.  Delta  Wing  Centerline 


Three-dimensional  effects  along  the  wing  centerline  of  delta  wings  at  angles  of 
attack  were  estimated  using  a  method  presented  in  References  2  and  4.  This  method 
is  based  on  numerical  solutions  by  the  method  of  Kennet  (Reference  12)  and  wedge 
theory.  The  method  of  Reference  12  is  valid  only  at  high  angles  of  attack  where  the 
shock  wave  is  detached  from  the  leading  edges,  while  wedge  theory  is  applicable  only 
at  low  angles  of  attack.  The  method  of  Reference  2  uses  the  results  of  Kennet  and 
wedge  theory  to  provide  a  means  of  making  estimates  at  intermediate  angles  of  attack. 

Xear  the  wing  centerline  it  was  found  that  the  edge  streamlines  could  be 
expressed  by: 
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A  correlation  for  n  reported  in  Reference  2  is  given  by: 
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is  presented  as  a  function  of  $**/£$  and  sweep  angle  A  in  Figure  9. 


In  estimating  three-dimensional  flou  effects  on  heating  it  is  again  necessary  to 
separate  the  influence  of  geometry  and  crossflow  pressure  gradients.  In  this  case  r 
can  vary  with  x  even  though  no  Ixxiy  curvature  exists.  This  variation  is  caused  by 
shock  wave  curvature,  and  cannot  be  easily  estimated.  However,  it  can  be  shown 
from  the  spanwise  momentum  equation  that  crossflow  pressure  gradients  will  cause 
streamline  curvature.  Hence,  it  is  concluded  that  when  the  streaml  nes  are  straight 
the  influence  of  crossflow  pressure  gradients  can  be  neglected  (df/dx  =  0).  This  con¬ 
dition  is  satisfied  when  n  is  either  zero  (two-dimensional  flow)  or  unity  (conical  flow). 
The  following  relationships  were  found  to  satisfy  these  two  conditions. 
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SECTION  VI 


NON1SOTH ERMAL  WALL  EFFECTS 


Nearly  all  methods  for  estimating  turbulent  heating  rates  are  applicable  only  to 
isothermal  surfaces.  However,  in  practice,  sizeable  wall  temperature  gradients  can 
exist  because  of  variations  in  local  heating  rates,  surface  emissivities,  and  the 
presence  of  internal  heat  sinks  (e.  g. ,  cryogenics).  In  evaluating  these  effects  it  is 
convenient  to  consider  separately  the  influence  of  variations  of  temperature  level  and 
that  of  temperature  gradients. 

Boundary  layer  growth  is  influenced  by  local  wall  ;emperature.  Consequently, 
the  boundary  layer  thickness,  and  therefore,  the  local  heat  transfer  coefficient,  is 
dependent  on  the  wall  temperature  at  all  upstream  locations.  The  influence  of  up¬ 
stream  variations  in  wall  temperature  are  reflected  through  the  Prpr  terms  in  the 
equivalent  distances  defined  by  Equations  (8)  and  (9).  However,  these  effects  on  heat 
transfer  rates  are  usually  much  smaller  than  the  effects  of  the  corresponding  thermal 
gradients. 

The  following  paragraphs  describe  a  method  presented  in  Reference  1  for  esti¬ 
mating  effects  of  thermal  gradients.  In  this  method  the  influence  of  thermal  gradients 
on  heat  transfer  rates  is  accounted  for  by  introducing  an  effective  thermal  potential 
term  $ ,  where: 


q  =  H  (i 
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(54) 


The  subscript  o  denotes  that  the  wall  enthalpy  is  evaluated  at  x  =  0.  The  method 
for  evaluating  $  is  essentially  a  modification  of  a  method  proposed  by  Seban  and 
described  in  Reference  13.  The  approximation  of  suggested  by  Seban  is: 


$ 

T, x=Xj 


(55) 


where 

x  =  x/x 


3  The  derivation  presented  in  Reference  13  is  based  on  wall  temperatures  instead  of 
enthalpies. 


In  spite  of  its  innocuous  appearance,  Equation  (55)  cannot  he  easily  integrated. 
Even  numerical  integration  is  difficult  because  of  the  singularity  occurring  when 
x  =  Xj.  An  alternate  formulation  is  presented  in  Reference  14  that  permits  an  easier 
numerical  solution  as  follows: 
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/  Pr  Mr  method 
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x’  is  the  length  of  the 
isothermal  surface  upstream 
of  the  heated  surface 


a)  General  Correlation 


b)  Typical  Comparison 

Figure  10:  MEASURED  AND  THEORETICAL  TURBULENT  HEATING  RATES  WITH 
LINEAR  WALL  TEMPERATURE  GRADIENT  ON  A  SHARP  CONE 


The  foregoing  analysis  is  restricted  to  flat  plate  flows.  The  influence  of  pressure 
gradients  and  three-dimensional  effects  are  included  bv  replacing  the  dimensionless 
streamwise  distance  "x  with  an  equivalent  distance  Sj.,  where 
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The  basis  for  this  definition  of  Sj  rests  primarily  on  a  similar  analysis  for  laminar 
flows  presented  in  Reference  1. 


The  set  of  equations  for  estimating  the  effects  of  wall  temperature  gradients  is 
given  in  finite  difference  form  by: 
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SECTION  VII 


NOSE  BLUNTNESS  EFFECTS 


The  flow  field  surrounding  a  vehicle  in  hypersonic  flight  is  dependent  on  the  nose 
geometry.  The  presence  of  a  blunt  nose,  for  example,  tends  to  increase  static  tem¬ 
perature  and  decrease  velocity  at  the  boundary  layer  edge.  This  effect,  which  can 
extend  many  diameters  downstream  of  the  nose,  can  cause  a  substantial  decrease 
(30  to  40  percent)  in  aerodynamic  heating  rates. 

Nose  bluntness  effects  are  dependent  on  the  vehicle  configuration,  Mach  number, 
Reynolds  number,  wall  cooling,  and  total  enthalpy  (real  gas  effects).  Two  limiting 
cases  are  immediately  recognized.  A  good  estimate  of  the  upper  bound  on  heating  can 
be  obtained  by  assuming  sharp  body  values  for  local  velocity  and  enthalpy.  Conversely, 
the  lower  limit  is  obtained  by  assuming  all  of  the  fluid  in  the  boundary'  layer  has  passed 
through  a  normal  shock  in  computing  local  flow  properties.  The  flow  conditions  at  the 
boundary  layer  edge  are  then  obtained  assuming  an  isentropic  expansion  from  the  stag¬ 
nation  to  the  local  pressure.  This  approach  is  restricted  to  equilibrium  or  frozen  flows. 

An  approximate  method  for  interpolating  between  the  upper  and  lower  heating  limits 
was  developed  during  the  present  study.  The  derivation  of  this  method  is  presented  in 
this  section,  and  data-theory  comparisons  are  presented  in  Volume  II  of  this  report. 


1.  SHOCK  ANGLE  EFFECTS  ON  STANTON  NUMBER 


The  first  step  in  the  present  analysis  was  to  determine  the  influence  of  the  shook 
angle  e  on  the  Stanton  number  Ngt.  Calculations  of  Ngt  were  made  using  the  PrVr 
method  for  edge  velocities  and  static  enthalpies  corresponding  to  several  shock  angles. 
Local  static  pressures  were  assumed  to  be  unaffected  by  bluntness.  The  local  flow 
properties  were  calculated  assuming  an  isentropic  expansion  behind  the  shock  wave. 
The  results,  shown  in  Figure  11,  indicate  that: 
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The  subscript  NS  denotes  that  the  evaluation  is  made  for  a  normal  shock  (<  =  0)  and 
Sh  denotes  that  the  evaluation  is  for  a  sharp  body  with  an  attached  shock. 
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Figure  11:  VARIATION  OF  HEAT  TRANSFER  COEFFICIENT  WITH 
SHOCK  ANGLE  ON  A  BLUNT  FLAT  PIATE 


2.  MASS  CONSERVATION  ANALYSIS 


Substituting  Equation  (67)  into  Equation  (65)  and  solving  for  sin  e  gives 
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Methods  for  obtaining  the  streamtube  width  A  are  presented  in  Section  V.  On  the 
hemisphere,  the  flow  is  assumed  to  be  axisymmetric.  At  the  tangent  point  TP  of  the 
hemisphere-wing  centerline 


ii  R  cos  0 
n  m 


The  angle  delta  is  the  angle  of  the  oody  surface  with  respect  to  the  free  stream 
flow.  Using  Equation  (69),  Equation  (68)  becomes 


3.  EVALUATION  OF  (6  -  6*) 

If  pressure  gradient  effects  are  neglected,  the  thickness  parameter  (6  -6*)  is 
easily  related  to  the  momentum  thickness  0  .  Following  the  practice  of  Beckwith  and 
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Gallagher  (Reference  14),  it  is  assumed  that  the  velocity  and  enthalpy  profiles  can 
be  related  to  the  corresponding  incompressible  profiles  by  the  following  transformation: 


Y  = 


(P/P  )dy 
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6  =  Y  t  (71) 

y=6 

The  velocity  ratio  u/ug  is  then  considered  to  be  a  function  only  of  Y/6,  and  indepen¬ 
dent  of  wall  cooling  and  Mach  number.  Thus: 


If  the  influence  of  pressure  gradients  is  neglected,  then  u/ue  «  ( Y/6 ) 1 ,  which  gives: 
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The  momentum  thickness  6  can  be  obtained  either  by  using  the  prpr  program 
described  in  Volume  III  or  from  Equation  (18). 


4.  EVALUATION  OF  BLUNT  BODY  CONSTANT  CR 

Using  Equations  (70)  and  (73),  Equation  (62)  becomes: 


(7. 


em  is  assumed  to  correspond  to  the  equivalent  attached  shock  angle.  Bluntness 
affects  p  .  u„,  8 ,  and  A.  However,  including  this  effect  would  require  a  tedious 

iterative  treatment.  Excellent  agreement  with  data  was  obtained,  however,  when 
these  parameters  were  computed  for  an  attached  shock  along  with  the  approximation 
that  CR  =  7. 
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SECTION'  VIII 


APPLICATIONS 


Equations  and  correlations  defining  the  PrMr  method  for  computing  turbulent 
heating  rates  are  presented  in  preceding  sections  and  in  Appendix  A  of  this  volume. 

In  this  section,  applications  of  this  method  in  estimating  turbulent  heating  on  a  blunt 
delta  wing  configuration  during  orbital  reentry  are  described.  Three  gas  models  are 
used  in  order  to  illustrate  real-gas  effects.  The  application  of  the  Prpr  method  in 
extrapolating  experimental  heat  transfer  data  obtained  from  ground  facilities  to  flight 
conditions  is  also  discussed. 


1.  GEOMETRY  AND  TRAJECTORY 

The  configuration  selected  for  this  analysis  is  a  delta  wing  with  a  hemispherical 
nose  cap  and  cylindrical  leading  edges.  The  dimensions  are  given  in  the  sketch  shown 
below. 


The  flight  path  selected  is  20,  000  feet  below  the  equilibrium  glide  trajectory  for  a 
W/SCl  of  150.  This  path  represents  a  severe  reentry  heating  trajectory.  The 
15-degree  angle  of  attack  approximates  the  condition  for  maximum  lift  to  drag  ratio 
for  this  geometry. 


2.  FLOW  FIELD  AND  GAS  PROPERTIES 

The  following  information  regarding  flow  field  and  gas  properties  is  require!  in 
order  to  estimate  turbulent  heat  transfer  coefficients  using  the  prpr  method: 

1.  surface  pressure  and  streamwise  pressure  gradient 

2.  velocity,  ug,  and  enthalpy,  ie,  at  the  boundary  layer  edge 

3.  temperature,  Tw,  or  enthalpy,  iw,  at  the  wall 


4.  streamline  divergence  parameters  r  and  f 

5.  diffusion  parameter  £ 

6.  viscosity,  p ,  partial  Prandtl  number,  <xr,  and  compressibility  factor- 
temperature  product,  ZT. 

The  flow  field  parameters  (items  1  through  4)  must  be  specified  at  all  points 
along  a  streamline.  The  gas  properties  (items  5  and  6)  are  usually  specified  as  a 
function  of  pressure  and  enthalpy. 

a.  Surface  Pressures 


Pressures  along  the  delta-wing  lower  surface  are  obtained  using  methods  pre¬ 
sented  in  Reference  2.  The  expression  for  sharp  delta  wing  pressure  coefficients  is: 


Cp  =  2  (sin  a) 


sin  (a  +  £  ) 
cos  £ 


The  shock  standoff  angle  £  is  obtained  using: 


5  =  70  - 


p„  \  ( tan  a 


i .  566 


,tan  A , 


(degrees) 


(75) 


(76) 


Equation  (76)  is  an  empirical  fit  to  numerical  flow  field  solutions  obtained  on  a  sharp 
delta  wing  using  the  method  of  Reference  12.  Results  from  this  method  are  in  good 
agreement  with  those  from  Reference  4. 

The  influence  of  nose  bluntness  is  estimated  using  the  following  equation  suggested 
by  van  Hise  (Reference  15). 
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P 


-  C, 


Sharp 


.  12  CD,n  .9 
7  (S/D)  '  2 

A  ACO 


(77) 


The  drag  coefficient  for  the  hemisphere  (Cp  n)  is  assumed  to  be  0.8,  which  is 
obtained  by  integrating  the  hemisphere  pressures  given  in  Reference  4. 


Pressures  on  the  leading  edge  are  computed  assuming  the  shock  wave  is  parallel 
to  the  surface.  For  this  case,  the  stagnation  line  pressure  is  equal  to  the  total  pres¬ 
sure  based  on  the  flow  component  normal  to  the  leading  edge.  Except  very  near  the 
nose,  real-gas  effects  on  pressures  were  found  to  be  small  at  the  flight  conditions 
considered,  and  are  neglected  in  this  analysis. 
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b.  Edge  Velocities 


The  velocities  at  the  boundary  laver  edge  for  the  sharp  delta  wing  case  are  com¬ 
puted  using  Equation  (78)  which  i*.  an  empirical  expression  presented  in  Reference  4. 


u 

e 


2 

a 

5600 


(78) 


Where  «  is  in  degrees. 

The  velocities  for  the  normal  shock  calculations  are  computed  assuming  an  isen- 
tropic  expansion  from  the  stagnation  point.  Bluntness  effects  on  heat  transfer  are 
estimated  from  the  sharp  and  normal  shock  results  using  Equation  (74). 

Consistent  with  the  parallel  shock  assumption,  stagnation  line  velocities  are 
obtained  by: 


-  =  cos  A  „ 

u„  eff 


(79) 


For  an  unyawed  condition,  the  effective  sweep  angle  (Aeff)  is  expressed  by: 


eff 


=  sin 


-1  /sin  A 


l 


sin  A  \ 
sin  a/ 


(80) 


c.  Wall  Temperatures 

The  wall  temperatures  are  assumed  to  be  constant  over  the  vehicle  at  any  given 
instant  of  time.  These  temperatures  are  given  as  a  function  of  velocity  and  altitude 
in  the  following  table. 


,oo,  ft/sec 

ALT,  ft 

T  * 

18,700 

187,  000 

3040 

16,000 

171,500 

2960 

14,000 

160,500 

2840 

12, 000 

149, 000 

2660 

10, 000 

137,500 

2400 

8,  000 

125, 000 

2100 

These  temperatures  are  I’epresentative  of  equilibrium  values  for  a  radiation  cooled 
surface. 
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d.  Streamline  Divergence  Parameters 

The  streamline  divergence  paramours  r  and  f  are  computed  using  Equations 
(52)  and  (53)  presented  in  Section  V. 

e.  Gas  M^del 


Most  of  the  heat  transfer  calculations  presented  in  this  section  are  for  air  in 
chemical  equilibrium  for  an  ideal  gas.  However,  computations  for  a  frozen  flow  are 
shown  at  the  enthalpy  for  peak  turbulent  heating  (V,*,  =  18,  700  fps).  Since  —juilibrium 
flow  implies  infinite  reaction  rates  and  frozen  flow  zero  reactions,  differences  in 
heating  rates  for  these  conditions  are  indicative  of  the  maximum  influence  of  reaction 
rates  on  turbulent  heating. 

1)  Chemical  Equilibrium 

The  turbulent  heat  transfer  estimates  for  equilibrium  flow's  were  obtained  using 
the  prpr  program  described  in  Volume  III.  The  diffusion  parameter  £.  is  computed 
from  Equation  (21),  and  the  stagnation  reference  viscosity  from  Equation  (19). 
Other  thermodynamic  and  transport  properties  are  essentially  the  same  as  those  pre¬ 
sented  in  Reference  16. 

2)  Frozen  Flow 

For  this  case,  the  gas  is  assumed  to  be  frozen  at  the  stagnation  point  composition. 
This  composition  and  the  corresponding  thermodynamic  and  transport  properties  are 
computed  using  the  equations  presented  in  Reference  16.  The  stagnation  point  gas 
composition  was  obtained  from  Reference  17.  The  flow  properties  on  the  wing  sur¬ 
face  are  based  on  an  isentropic  expansion. 

For  the  range  of  temperatures  considered,  the  vibrational  energy  level  of  the 
diatomic  molecules  can  be  considered  to  be  proportional  to  the  translational  and 
rotational  energy  levels.  This  approximation  greatly  simplifies  computations  since 
the  specific  heats  Cp  and  cv  are  then  constant  with  temperature. 

By  definition  of  frozen  flow,  no  chemical  recombination  of  atoms  occurs.  Conse¬ 
quently,  the  energy  absorbed  in  dissociation  must  be  considered  unavailable  in  com¬ 
puting  the  total  enthalpy  of  the  flow’.  For  the  present  case,  the  total  enthalpy  of  the 
flow  at  the  boundary  layer  edge  is  reduced  by  dissociation  to  about  47  pert  ent  of  the 
free-stream  value  at  Vw  =  18,700  fps.  The  resulting  decrease  in  adiabatic  w’all 
enthalpy  causes  a  substantial  reduction  in  heating  rates.  The  frozen  flow  adiabatic 
wall  enthalpy  is  given  by: 

i  -  rl  ■+  (1  -  r)i  =  .  47  rl  -*  (1  -  r)i  (81) 

aw  e  e  00  e 


If, 


A  second  frozen-flow  heat-transfer  calculation  was  made  that  differs  from  that 
described  previously  only  in  that  the  adiabatic  wall  enthalpy  is  based  on  the  total 
enthalpy  of  the  free-stream  flow  as  follows: 

i  =  rl+  (1  -  r)i  (82) 

The  second  approach  is  presented  to  furnish  an  estimate  of  the  influence  of  sudden 
recombination  of  atoms  or.  heating.  A  noneatalytic  wall  is  assumed  in  all  cases. 

3)  Ideal  Gas 

Since  the  gas  composition  is  assumed  to  be  frozen  in  the  free-stream  state,  the 
ideal  gas  case  is  a  special  type  of  frozen  flow.  Viscosities  are  computed  using 
Sutherland's  law.  The  usual  thermodynamic  parameters  for  low  temperature  air 
given  below  were  assumed. 


y  =  1.4 

9  9 

Cp  =  6006  ft  /sec  -°R 
R  =  1716  ft2/sec2-°R 
Z  =  1.0 


3.  BLUNTNESS  AND  REAL  GAS  EFFECTS 

Bluntness  effects  on  the  Stanton  number  along  the  wing  centerline  are  illustrated 
in  Figure  12.  At  the  highest  velocity  considered,  the  influence  of  bluntness  is  seen  to 
extend  only  about  2  diameters  downstream  of  the  nose  cap.  At  lower  velocities  and 
altitudes,  the  boundary  layer  displacement  thickness  is  smaller;  consequently,  blunt¬ 
ness  effects  extend  much  further  downstream.  The  sharp  body  estimates  in  computing 
the  blunt  body  heating  rates  are  based  on  local  velocities  and  enthalpies  for  a  sharp 
wii.g,  but  are  adjusted  to  account  for  bluntness  effects  on  local  pressure. 

Real  gas  effects  are  illustrated  in  Figure  13  which  shows  Stanton  numbers  for 
equilibrium,  frozen,  and  ideal  gas  flows,  all  based  on  normal-shock  theory  (see 
Section  VI).  The  frozen  flow  estimates  are  seen  to  be  substantially  lower  than  the 
equilibrium  values,  even  when  the  adiabatic  wall  enthalpy  is  based  on  the  free-stream 
total  enthalpy. 

The  comparisons  shown  in  Figure  13  are  unrealistic  in  that  the  frozen  flow  com¬ 
position  applies  only  to  the  flow  originating  at  the  stagnation  point.  A  more  meaningful 
comparison  is  given  in  Figure  14,  which  shows  the  blunt  body  estimates  obtained  using 
Equation  (74),  The  frozen  flow  composition  is  assumed  in  computing  the  normal-shock 


Turbulent  flow 


S/D 

Figure  12:  NOSE  BLUNTNESS  EFFECT  ON  DELTA  WING 
CENTERLINE  STANTON  NUMBER 


Figure  14:  REAL-GAS  EFFECT  ON  DELTA-WING  CENTERLINE 
STANTON  NUMBER  USING  BLUNT -BODY  METHOD 


values,  and  equilibrium  flow  is  used  in  computing  the  sharp  wing  heating  rates.  It  is 
seen  that  for  this  case  reaction  rates  influence  heating  only  near  the  nose  cap.  The 
ideal-gas  Stanton  numbers  are  about  40  percent  lower  than  the  corresponding  equilib¬ 
rium  values. 


4.  EXTRAPOLATION  TO  FLIGHT 

At  this  time  no  test  facility  except  hypersonic  gun  ranges,  is  capable  of  duplicating 
hypersonic  flight  conditions  for  speeds  greater  than  about  10,  000  feet  per  second. 
Except  for  shock  tubes,  the  total  enthalpy  is  less  than  high-speed  flight.  Shock  tubes 
are  capable  of  simulating  enthalpy  but  are  limited  to  low  Mach  numbers.  For  this 
reason,  empirical  methods  alone  cannot  be  relied  upon  to  provide  accurate  heat  trans¬ 
fer  predictions  for  hypersonic  flight.  However,  analytical  methods  alone  are  not  ade¬ 
quate  for  predicting  heating  on  realistic  aerodynamic  configurations  unless  supported 
by  test  data  for  that  specific  configuration.  It  therefore  becomes  necessary  to  inter¬ 
pret  test  data  obtained  from  ground  facilities  in  such  a  way  that  these  data  can  be 
reflected  in  predictions  for  flight.  One  approach  for  making  this  interpretation  is  to 
extrapolate  these  data  to  flight  conditions  using  analytical  methods.  In  performing 
this  extrapolation,  it  is  convenient,  whenever  possible,  to  normalize  Stanton  numbers 
with  respect  to  some  reference  condition  such  that  the  normalized  Stanton  numbers 
are  not  greatly  dependent  on  flow  conditions.  The  reference  condition  selected  for 
this  analysis  is  the  stagnation  line  of  an  infinite  cylinder  of  one  foot  diameter  and  a 
sweep  angle  of  60  degrees.  The  reference  Stanton  numbers  for  the  flight  path  being 
considered  is  presented  in  Figure  15  for  both  equilibrium  and  ideal-gas  flows. 

The  Stanton  number  ratios  at  the  wing  centerline  are  shown  in  Figure  16.  Similar 
comparisons  for  the  leading  edge  stagnation  line  and  the  lower-surface  leading-edge 
tangent  lines  are  presented  in  Figure  17.  Extrapolation  factors  for  two  typical  wind- 
tunnel  test  conditions  are  shown  in  Figure  18.  The  good  simulation  of  Stanton  number 
ratios  for  the  stagnation  line  is  not  surprising  since  the  reference  Stanton  number  is 
also  for  the  stagnation  line  but  at  a  different  sweep  angle.  However,  significant  differ¬ 
ences  in  Stanton  numbei  ratios  are  noted  for  the  delta  wing  centerline,  particularly  for 
Arnold  center  tunnel  B.  These  deviations  arc  due  to  differences  in  local  Mach  number, 
enthalpy  level,  and  the  ratio  of  the  wall  enthalpy  to  the  total  enthalpy. 


Figure  15:  REAL-GAS  EFFECT  ON  REFERENCE  STANTON  NUMBER 


Figure  16:  REAL-GAS  EFFECT  ON  DELTA  WING  CENTERLINE  STANTON 
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Figure  17:  STANTON  NUMBER  RATIOS  ON  THE  CYLINDRICAL 
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SECTION  IX 


CONCLUDING  REMARKS 


The  prpr  method  for  estimating  turbulent  heating  rates  is  well  substantiated  by 
experimental  data  obtained  on  several  body  shapes  and  for  a  wide  range  of  test  condi¬ 
tions.  In  particular,  effects  of  local  Mach  number,  wall  cooling  ratio  (le/iw),  three- 
dimensional  flow,  and  streamwise  pressure  gradients  predicted  by  this  method  are  in 
excellent  agreement  with  experimental  trends.  Similarly,  the  method  for  estimating 
nose  bluntness  effects  described  in  Section  VII,  in  general,  furnishes  good  agreement 
with  available  heat  transfer  data. 

Methods  for  predicting  real-gas  effects  are  not  as  well  established.  The  combi¬ 
nation  of  high  pressures  and  temperatures  required  to  obtain  a  highly  dissociated 
turbulent  flow  is  difficult  to  achieve  in  ground  facilities;  consequently,  very  little  data 
are  available  for  turbulent  flows  with  significant  levels  of  dissociation  and  high  Mach 
numbers. 

Whenever  possible,  it  is  recommended  that  the  prpr  computer  program  described 
in  Volume  III  be  used  for  calculating  convective  heat  transfer  rates.  However,  if  com¬ 
putations  must  be  made  manually,  the  simplified  equations  presented  in  Appendix  B  are 
recommended.  In  most  cases,  the  simplified  equations  provide  heat  transfer  estimates 
within  a  few  percent  of  the  complete  equations.  Handbook  methods  presented  in  Refer¬ 
ence  5  will  also  provide  heating  estimates  that  are  in  good  agreement  with  the  pr  pr 
computer  program  except  for  delta  wing  surfaces.  Significant  discrepancies  in  heating 
estimates  on  sharp  delta  wings  can  occur  because  of  differences  in  methods  of  deter¬ 
mining  three-dimensional  parameters. 

The  purpose  of  developing  the  turbulent  nonsimilar  method  was  to  furnish  a  basis 
for  evaluating  the  profile  parameters  appearing  in  the  equations  and  to  substantiate 
the  method  presented  in  Section  VI  for  estimating  effects  of  wall  temperature  gradients. 
Unfortunately,  it  was  impossible  to  complete  these  studies  during  the  present  investi¬ 
gation.  Numerical  instabilities  often  developed  when  pressure  gradients  were  imposed 
on  the  flow.  Although  measures  were  usually  found  for  eliminating  these  instabilities, 
they  resulted  in  substantial  increases  in  computing  time. 

In  spite  of  these  difficulties,  the  turbulent  nonsimilar  method  is  believed  to  repre¬ 
sent  a  new  and  promising  approach  in  treating  turbulent  flows.  This  approach  offers 
several  advantages  over  previous  methods.  Empirical  correlations  are  used  only  in 
defining  Reynolds  stresses,  and  no  assumptions  regarding  velocity  profiles  or  shear 
distributions  are  required.  Since  the  flow'  equations  are  solved  in  partial  differential 
form,  no  coordinate  transformations  are  required.  Such  transformations  usually 
impose  restrictions  on  boundary  conditions  and  flow  similarity.  Skin-friction  coeffi¬ 
cients  and  heat-transfer  rates  obtained  using  this  method  are  in  good  agreement  with 
experimental  data  over  a  wide  range  of  test  conditions. 


APPENDIX  A 


DERIVATION  OF  THE  p  g  METHOD 

rr  r 

The  derivation  of  the  p  u  equations  presented  in  this  appendix  is  based  on  a 

r  r 

solution  of  the  boundary  layer  energy  integral  equation.  A  similar  derivation  based 
on  a  solution  of  the  momentum  integral  equation  is  presented  in  Appendix  B  of 
R'i  "erenee  2.  The  present  derivation  is  given  in  the  following  parts: 

1.  Derivation  of  a  general  form  of  the  boundary  layer  energy  integral  equation. 

2.  Transformation  and  solution  of  the  energy  integral  equation. 

3.  Solution  of  the  momentum  integral  equation. 

4.  Evaluation  of  laminar  boundary  layer  parameters. 

5.  Combined  laminar  and  turbulent  method. 

A  simplified  version  for  making  hand  calculations  is  presented  in  Appendix  B  and 
the  I3M  7094  digital  computer  program  using  the  equations  presented  in  this  appendix 
is  described  in  Volume  m. 

The  correlations  given  in  part  4  of  this  appendix  are  the  same  as  those  presented 
in  Reference  1,  but  differ  from  those  of  Reference  2.  However,  differences  in  heating 
rates  computed  using  the  present  correlations  and  those  of  Reference  2  are  small. 


1.  DERIVATION  OF  THE  ENERGY  INTEGRAL  EQUATION 

This  derivation  is  restricted  to  the  vieinity  of  a  plane  of  symmetry  as  well  as  by 
the  usual  boundary  layer  assumptions.  The  mass  and  energy  conservation  equations 
for  boundary  layer  flows  are  given  by  Equations  (A-l)  and  (A-2),  respectively: 
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The  corr  .'spending  control  volume  is  shown  in  the  following  sketch: 


z,w 


The  length  elements  in  x  and  z  are  unity.  However,  the  length  element  in  y  is 
determined  by  the  function  Ay  =  g(x),  which  remains  arbitrary  (subject  to  the  restric¬ 
tion  that  dg/dx  remains  finite).  Later  it  will  be  seen  that  in  most  cases  the  most 
convenient  choice  of  g  is  dependent  on  the  geometry  of  the  body  under  consideration. 
The  surface  y  =  0  is  by  definition  a  line  of  symmetry,  hence  v.  but  not  necessarily 
3v/3y  ,  is  zero  when  y  =  0.  Also,  it  should  be  noted  that  for  turbulent  flows  the  flow 
parameters  appearing  in  Equations  (A-l)  and  (A-2)  are  averaged  with  respect  to  time. 
For  turbulent  flows,  t  represents  the  effective  shear  stress  including  the  Reynolds 
stress,  and  q  represents  a  corresponding  effective  heat  flux. 

The  velocity  component  normal  to  the  body,  w  ,  is  found  by  integrating 
Equation  (A-l)  with  respect  to  z: 


Using  Equation  (A-3),  and  assuming  that  w  =  0,  Equation  (A-2)  can  be  written: 
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I..  *  Jer  to  obtain  the  boundary  layer  momentum  integral.  Equation  (A-4)  is  integrated 
from  the  wall  to  some  arbitrary  location  (z  =  h)  outside  of  the  boundary  layer. 
Neglecting  external  vorticity,  q  =  r  ^  =  0,  the  momentum  integral  becomes: 
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Integrating  the  second  term  by  parts  and  rearranging,  Equation  (A-5)  can  be 
expressed: 

V-0  P"*(Ie  -  »  dz  +  i  jf  pV(Ie  ‘  »  dZ  (A'6) 

Since  this  analysis  is  restricted  to  lines  of  symmetry,  it  can  be  easily  shown  that: 
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In  the  present  analysis  it  is  further  assumed  that: 
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The  influence  of  wall  temperature  gradients  on  heal  transfer  is  treated  separately  in 
Section  Vi. 


Thus.  Equation  (A-G)  can  now  be*  expressed  in  terms  of  the  heat  transfer  coefficient, 
H,  as  follows: 


(A-7) 


Now,  introducing  the  boundary  layer  thickness  parameters,  energy  thickness: 


and  the  crossflow  energy  thickness  ratio: 


Equation  (A-7)  now  reduces  to: 
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Note  that  e  and  9ve/9y  are  determined  by  the  definition  of  g,  which  is  still 
arbitrary.  Since  results  obtained  from  Equation  (A-9)  are  independent  of  this  defini¬ 
tion,  the  selection  is  made  solely  on  the  basis  of  convenience  in  evaluating  the 
appropriate  flow  parameters. 


It  is  seen  from  Equation  (A-8)  that  an  exact  evaluation  of  e  requires  a  solution  of 
both  the  streamwise  and  crossflow  velocity  profiles.  Such  solutions  are  available  only 
for  laminar  flow  over  simple  shapes.  It  is  therefore  desirable  to  define  g  in  such  a 
way  that  e  and  9\’e/9y  reflect  only  the  influence  of  crossflow  pressure  gradients. 

For  this  case  g  reflects  only  the  streamline  divergence  due  to  body  geometry.  In  the 
absence  of  ^  more  exact  approach,  g  can  be  assumed  to  be  proportional  to  the  body 
radius  of  curvature,  r,  normal  to  the  streamline.  In  order  to  be  consistent  with  the 


most  common  symbology,  r  is  used  in  place  of  g  in  the  following  analysis,  Also,  the 
streamline  divergence  due  to  crossflow  pressure  gradients  is  denoted  by  f,  where  f 
is  defined  by: 
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Using  Equation  (A-10),  Equation  (A-9)  can  now  be  expressed  as: 
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2.  TRANSFORMATION  AND  SOLUTION  OF  THE  ENERGY  INTEGRAL  EQUATION 
a.  Transformation 

In  order  to  obtain  a  more  useful  form  of  the  energy  equation,  a  modified 
Stewartson  transformation  suggested  by  Mager  (Reference  3)  is  adopted  in  which: 
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The  stagnation  values  of  density  and  viscosity,  p0  and  p0,  are  required  to  be 
constant,  and  F  is  an  unspecified  function  of  X  only.  With  these  definitions,  the 
energy  thickness  and  heat  transfer  coefficient  in  the  transformed  coordinate  system 
are,  respectively: 
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c  ,  r,  and  f  are  unchanged  by  the  transformation.  The  transformed  energy  equation 
now  becomes: 
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rf 

b.  Solution  of  the  Transformed  Integral  Equation 


In  Appendix  B  of  Reference  2,  the  transformed  momentum  integral  equation  was 
solved  by  assuming  a  transformed  Blasius  shear  law  given  below: 
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The  exponent  m  is  unity  for  laminar  flows,  and  approximately  4  for  turbulent 
boundary  layers.  However,  by  leaving  m  unspecified  the  following  analysis  is  valid 
for  both  types  of  flow. 

Assuming  unit  Prandtl  number  and,  neglecting  the  effects  of  streamwise  pressure 
gradients,  the  enthalpy  can  be  related  to  velocity  by  the  well-known  Crocoo  energy 
relationship  given  by; 
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In  the  untransformed  physical  plane,  Equation  (A- 19)  is: 


where  Cx  is  approximately  0.  332  for  laminar  flow,  and  is  a  constant  to  be  specified 
later  for  turbulent  flow.  Using  Equation  (A-22),  H  can  now  be  expressed  in  the  more 
familiar  form  shown  below: 
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where  xeq  is  an  equivalent  flat  plate  length  parameter  defined  by: 
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Equations  (A-23)  and  (A-24)  are  equivalent  to  the  result  reported  in  Appendix  B  of 
Reference  2  with  one  exception.  The  exponent  ?  in  Equation  (A-24)  appears  in  place 
of  the  crossflow  momentum  thickness  ratio,  E  ,  in  Reference  2,  where: 


(A-25) 


Tf  the  Crocco  energy  relationship  is  assumed,  it  is  easily  seen  that  e  and  E  are 
equal.  Considering  that  approximate  methods  must  be  used  in  evaluating  either  e  or 
E  ,  particularly  for  turbulent  flow,  it  appears  reasonable  to  have  the  same  approxima¬ 
tions  for  both  parameters.  Accordingly,  e  is  replaced  by  E  in  the  following  analysis 
is  order  to  make  the  nomenclature  consistent  with  that  of  previous  publications  (e.  g.  . 
References  1,2,  and  3). 


3.  SOLUTION  Or  THE  MOMENTUM  INTEGRAL  EQUATION 

The  purpose  of  this  derivation  is  to  obtain  expressions  for  momentum  thickness 
and  skin  friction  consistent  with  the  heat  transfer  expression  given  by  Equations 
(A-23)  and  (A-24). 


05 


The  boundary  layer  momentum  integral  is  given  in  Appendix  B  of  Reference  2  by: 
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where  0  is  the  boundary  layer  momentum  thickness. 

Equation  (A-26)  can  be  rearranged  into  a  simple  form  similar  to  the  energy 
integral  expression  given  by  Equation  (A-ll),  giving: 
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where  A  =  2  +  ^*/0  . 

Applying  the  Magei.  transformation  defined  by  Equation  (A-12),  Equation  (A-27) 
becomes: 
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By  replacing  the  left  side  of  Equation  (A-28)  with  the  transformed  Blasius  shear  law, 
Equation  (A-15),  an  expression  for  ®  is  obtained  in  the  same  way  as  was  Q  in  part  2 
of  this  appendix.  The  result  is: 


®  = 


X. 


/ 

I/O 


m+1 

nw-1  „  „A-1  A  1/m  £  m 

-  C  F  p  U  u  (rP)  dX 

mm  ro  e  o  '  ' 


m 

m+1 


Fp  U  rf1 
o  e 


or, 


(A-29) 


In  untransformed  physical  coordinates,  the  equivalent  expression  for  8  is: 
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The  momentum  thickness  is  related  to  skin  friction  by  the  modified  Blasius  shear 
law  given  by  Equation  (A-15).  The  untransformed  equivalent  of  Equation  (A-15)  is 
given  by: 
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Substituting  Equation  (A-30)  into  (A-31)  gives: 
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Assuming  that  Cm  is  independent  of  x ,  it  can  be  expressed  in  a  form  similar  to  that 
for  Cq  in  Equation  (A-22)  by  the  following: 
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the  profile  parameter  Cx  is  the  flat  plate  value  for 
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and  C 


and  P  reflects  the  influence  of  streamwisc  pressure  gradients. 


Using  Equation  (A-33),  Equation  (A-32)  can  be  reduced  to  a  form  similar  to  that 
for  H  given  by  Equation  (A-23).  The  corresponding  skin  friction  expression  is: 
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where: 
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and  the  skin-friction  equivalent  distance  is  defined  by: 
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Using  these  definitions  for  Nn  „  c  and  ,  Equation  (A-30)  simplifies  to: 
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4.  EVALUATION  OF  L i  4INAR  BOUNDARY  LAYER  PARAMETERS 

Exact  solutions  c  the  similarity  form  of  the  laminar  boundary  layer  equations 
were  used  to  evaluate  the  parameters  appearing  in  Equations  (A-23),  (A-24)  and 
(A-32).  This  was  done  in  an  orderly  manner,  beginning  with  two-dimensional 
constant-property  constant-pressure  flow,  and  progressing  to  the  most  complex  con¬ 
ditions  for  which  exact  solutions  are  available.  The  evaluations  determined  from  the 
simpler  cases  were  retained  or  amplified  in  analyzing  the  more  complex  cases.  Thus, 
the  constant  Cx  for  laminar  flow  is  always  taken  to  be  0.  33206.  the  value  given  by 
Howarth  in  Reference  18  for  incompressible  flat  plate  flow.  The  effects  of  pressure 
gradients,  wall  cooling,  etc.  ,  are  accounted  for  in  other  terms  of  Equations  (A-23) 
and  (A-32). 
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In  some  cases,  alternate  definitions  were  possible.  For  example,  the  authors  of 
Rcfeiences  19  and  20  incorporated  (in  effect)  pressure  gradients  into  the  term  prpv 
appearing  in  Equation  (A-23),  while  in  the  present  formulation  such  effects  appear  in 
the  heat  transfer  equivalent  distance  x  .  The  latter  definition  is  to  be  preferred  as 
the  former  cannot  be  made  consistent  with  the  results  of  Reference  12,  which  presents 
solutions  for  various  pressure  gradients,  but  with  pp  held  constant.  The  definitions 
used  here  were  adopted  only  after  an  examination  of  several  possible  alternatives. 

The  criteria  for  selection  were  consistency  between  the  results  of  the  various  special 
cases,  conformity  with  physical  considerations,  accuracy,  simplicity,  and  freedom 
from  interdependencies. 

a.  General  Considerations 


As  a  matter  of  physical  consistency,  it  is  required  that  if  the  fluid  properties  p 
and  p  are  constant  through  the  boundary  layer  then  the  reference  values  of  the  fluid 
properties  be  equal  to  those  constant  values.  This  principle  is  extended  to  constant 
products  as  well,  i.  e.  ,  it  is  required  that  when  in  a  given  numerical  calculation 
(e.  g.  ,  References  12  and  21)  the  product  of  density  and  viscosity  is  held  constant  at 
some  base  value  (usually  the  wall),  then  the  reference  density-viscosity  product, 
prpr  ,  must  also  be  equal  to  that  base  value.  The  functions  Fpr  and  £  are  equal  to 
1.  0  when  a  and  Nj^  are  equal  to  1.0,  and  £  =  1.  0  for  ideal  gases.  Also,  in  flat 
plate  flow  the  equivalent  distance  is  equal  to  the  physical  distance  from  the  leading 
edge. 

b.  Two-Dimensional  Flat  Plate  Flow 

The  special  case  of  two-dimensional  flat  plate  flow  is  examined  first  since  the 
effect  of  fluid  property  variations  within  the  boundary  layer  can  be  examined  without 
the  additional  complexity  of  streamwise  variations.  For  the  case  of  constant  fluid 
properties,  the  solutions  of  Howarth  (Reference  18)  show  that  m  =  1  and  Cx  =  .  332, 
so  that  Equation  (A-23)  becomes: 
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where  Equation  (A-39)  follows  from  the  principles  stated  in  the  preceding  paragraph. 
For  this  special  case,  the  only  undetermined  quantities  are  the  Reynolds  analog}' 
factors  £  and  FPr  1  Note  that  the  reference  stagnation  viscosity,  p0  .  no  longer 
appears. 


(!9 


c. 


Reynolds  Analogy  Factors 


Following  the  practice  of  Reference  19,  for  example,  the  Prandtl  number  effect 
is  correlated  in  terms  of  a,  the  partial  Prandtl  number  for  translation,  rotation,  and 
vibration.  The  Prandtl  number  effect  on  Reynolds  analogy  in  flat  plate  flow,  usually 
given  as  Fpr  =  (a)2^3  for  constant  c,  is  somewhat  better  represented  by  cr*  as 
may  be  seen  in  Figure  19. 

For  variable  Prandtl  number  there  is  an  uncertainty  as  to  which  value  should  be 
used  in  correlating  its  effect.  All  solutions  in  the  literature  for  which  the  Prandtl 
number  is  variable  also  involve  variable  pp,  so  that  prpr  is  not  necessarily  equal 
to  pepe.  For  such  cases,  it  was  found  that  the  Prandtl  number  should  be  evaluated 
at  the  enthalpy  and  pressure  corresponding  to  prpr.  This  value  of  the  Prandtl 
number  is  hereafter  denoted  as  crr  . 

When  the  Prandtl  number  effect  is  correlated  in  terms  of  the  partial  Prandtl 
number,  the  effect  of  energy  transport  by  diffusion  must  be  treated  separately.  The 
first  exact  calculation  of  this  effect  is  reported  in  Reference  19,  wherein  the 
expression: 
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was  found  to  agree  well  with  exact  solutions  for  Nj^  =1.4  in  stagnation  point  flow. 
In  high  Mach  number  flows,  however,  Equation  (A-40)  may  predict  a  significant 
diffusion  effect  under  conditions  for  which  no  dissociation  actually  exists,  since  the 
temperatures  within  the  boundary  layer  may  be  well  below  the  stagnation  value.  To 
avoid  this  inconsistency,  Equation  (A-40)  was  modified  to  operate  on  the  local 
reference  enthalpy,  corresponding  to  prpr,  rather  than  the  stagnation  value.  The 
modified  expression  is  given  by: 
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Calculations  made  by  Cohen  (Reference  22)  of  Lewis  number  effects  for  non¬ 
stagnation  boundary  layers  are  shown  in  Figure  20,  along  with  computations  from 
Equation  (A-41). 

d.  Reference  Density-Viscosity  Product 

The  reference  density-viscosity  product  was  first  evaluated  for  zero  Mach  number 
with  various  degrees  of  wall  cooling  using  the  solutions  of  References  20,  23,  and  24 
and  some  unpublished  solutions  by  Halvorson  and  Cassmeyer  of  The  Boeing  Company, 
as  shown  in  Figure  21. 
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Figure  19:  LAMINAR  PRANDTL  NUMBER  EFFECT  ON 
REYNOLDS  ANALOGY  FACTOR 
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Figure  21:  CORRELATION  OF  REFERENCE  DENSITY-VISCOSITY  PRODUCT 


For  edge  Mach  numbers  greater  than  zero,  it  was  found  that  the  reference 
density-viscosity  product  prpr  can  be  represented  as  a  function  only  of  PeMe  , 
pwpw  and  pgi p s'  where  the  latter  is  the  density- viscosity  product  evaluated  at 
stagnation  enthalpy  and  the  local  pressure.  The  correlation  obtained  is  g’ven  below: 


where 


p  u  =  p  u 
Kmr  w*w 


PWPW 


7  , 

8 


1.2 


.2  + 


PWPW 


1_ 

10 


(A-42) 


and 


(P  P  )  rr  = 

e  e  eff 


K* 


13  3  /  „* 

16+16\  PwPwy 


pepe 


K*  = 


PS,PS' 

pepe 


1.005 


,1 

14 


005  + 


^PS,MS' 

pepe, 


Equation  (A-42)  is  plotted  in  Figure  2  (Section  III.  1)  using  the  approximation  that  the 
exponent  K*  f(pepe)/(PwPw)]  is  unity.  The  contribution  of  this  term  is  small,  and 
Figure  2  can  be  used  to  obtain  prpr  without  introducing  significant  errors.  Sub¬ 
sequent  investigations  have  shown  that  PrPr  is  independent  of  pressure  gradients,  as 
demonstrated  by  comparisons  with  exact  solutions  shown  in  Figure  22. 

e.  Pressure  Gradient  Effects 


1)  Evaluation  of  and  E 

Referring  to  Equation  (A-23)  and  recalling  the  earlier  comment  that  Cx ,  m, 

Fpr  ,  and  «£  are  by  definition  taken  as  the  flat  plate  values,  it  is  seen  that  all  pres¬ 
sure  gradient  effects  are  reflected  in  PrPr  and  xeq.  These  effects  can  be  evaluated 
for  similar  flows  from  the  solutions  published  which  consider  streamwise  pressure 
gradients,  (e.  g.  ,  References  12  and  20)  and  those  which  consider  crossflow  pressure 
gradients  (e.  g.  ,  References  21  and  22) 
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Figure  22:  REFERENCE  DENSITY-VISCOSITY  PRODUCT  FOR  FLOWS  WITH  PRESSURE  GRADIENT 


Beginning  with  the  simplest  possible  case,  two-dimensional  flow  of  an  ideal  gas 
writh  unit  Prandtl  number,  and  the  viscosity  proportional  to  temperature,  the  equiva¬ 
lent  distance  effects  can  be  isolated.  Since  pp  is  always  equal  to  pepe  ,  then  prPr 
is  also  equal  to  p^pg.  (Note  that  pp  is  not  necessarily  constant  through  the  flow 
field,  but  varies  with  the  local  boundary  layer  edge  pressure. )  With  these  values 
incorporated,  the  equivalent  distance  expression  in  Equation  (A-24)  reduces  to: 
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In  Equation  (A-43)  the  term  pepeug  reflects  the  effects  of  upstream  variations,  while 
JL  accounts  for  local  streamwise  pressure  gradient  effects  on  the  boundary  layer 
profiles. 

A  correlation  has  been  f~und  for  JL  which  may  be  written: 
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where  0  is  the  dimensionless  pressure  gradient  parameter  similar  to  that  defined  by 
the  authors  of  Reference  11.  The  parameters  0,  ,  and  are  defined  by 

Equations  (A-50)  to  (A-53).  Subsequent  investigations  of  exact  solutions  for  non-unit 
Prandtl  number  and  nonlinear  viscosity  laws  have  shown  that  expressions  of  the  form 
of  Equation  (A-43)  are  valid  for  these  more  complex  conditions  as  well,  either  for 
two-dimensional  flows  with  streamwise  pressure  gradients,  or  for  yawed  cylinder 
flow.  The  expressions  finally  adopted  are: 
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and  a  generalization  of  Equation  (A-43): 
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exp  K  =  0  when  N  <  .  05  and  .  99  <  N  <  1.  01 , 

exp  K  =  194  exp  [  -f  N(N-l)  I  when  .  05  <  N  <  .  99, 

l  O  J 

exp  K  =  .  194  exp  J  (N-3 '  j  when  N  >  1.  01 ,  and  N  =  ^ 

The  subscripts  s  and  c  are  introduced  to  distinguish  between  streamwise  and 
crossflow  pressure  gradients;  it  should  be  noted  that  JL  is  concerned  only  with 
streamwise  pressure  gradient  effects  and  only  with  crossflow  effects;  also,  note 
that  -  1.0  for  Ps  =  0. 


The  functions  F^.  (3 ,  F^  ,  and  Z  are  given  for  either  streamwise  or  crossflow 
pressure  gradients  by  the  following  expressions: 

(a)  Streamwise  Pressure  Gradients: 
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where 
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The  subscript  m  denotes  evaluations  at  local  pressure  and  a  mean  boundary  layer 
enthalpy  defined  by: 
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The  second  equality  in  Equation  (A-53)  follows  from  the  condition  of  constant  pressure 
across  the  boundary  layer. 

(b)  Crossflow  Pressure  Gradients: 
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The  subscript  SL  refers  to  the  stagnation  line  and  (ZT)m  is  evaluated  at  local 
pressure  and  a  mean  boundary  layer  enthalpy  defined  by: 
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With  a  minor  modification  of  xetj  it  is  easily  shown  that  the  definition  of  0S  , 
given  by  Equation  (A-51).  is  identical  to  the  correspond'  » leter  of  Reference  12. 

The  expressions  for  obtaining  Jl.  Equations  (A-48)  ana  (A- 50)  to  (A-54)  were 
developed  on  the  basis  of  providing  the  best  fit  to  the  exact  similar  solutions  shown  in 
Figure  23.  However,  the  obvious  similarity  of  Equations  (A-50)  to  (A-54)  to  the 
various  reference  enthalpies  appearing  in  the  literature  provides  some  analytical 
justification  for  these  correlations. 


The  equations  for  obtaining  J ^  furnished  a  basis  for  determing  E.  The  cross¬ 
flow  pressure  gradient  parameter  /?c  was  assumed  to  be  unity,  which  is  the  value  for 
an  unyavved  cylinder.  Comparisons  of  E  from  Equations  (A-49)  and  (A-55)  to  (A-58) 
for  swept  cylinder  flow  are  shown  in  Figure  24.  Similar  comparisons  with  yawed 
cone  solutions  are  presented  in  Figure  25.  The  term  (2Zc)exPK  was  developed  on 
the  basis  of  the  yawed  cone  comparisons.  Note  that  this  term  is  unity  for  the  cylinder 
case. 


2)  Evaluation  of  PL 

Correlations  of ’’exact"  laminar  solutions  show  that  the  influence  of  streamwise 
pressure  gradients  on  skin  friction,  which  is  reflected  in  Pj^,  can  be  easily  related 
to  the  corresponding  effect  on  Jj^.  The  relationship  is: 


P 


L 


(A-59) 


Comparisons  with  solutions  presented  in  References  11  and  22  are  shown  in  Figure 26. 
Equation  (A-36)  now  becomes: 


S  = 
cq 


(A-60) 
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Figure  24:  REAL-GAS  CORRELATION  PARAMETER  FOR  SWEPT  INFINITE  CYLINDERS 


Figure  26:  STREAMWISE  PRESSURE  GRADIENT  EFFECT  ON  SKIN  FRICTION 


5.  COMBINED  LAMINAR  AND  TURBULENT  METHOD4 

The  PrMr  computer  program  described  »n  Appendix  D  is  intended  for  making 
heat  transfer  and  skin  friction  predictions  for  both  laminar  and  turbulent  flows.  For 
computer  purposes,  it  is  convenient  tc  utilize  parameters  common  to  both  laminar  and 
turbulent  boundary  layers  whenever  possible.  A  problem  is  encountered  in  defining  a 
reference  Reynolds  number  since  the  equivalent  distances  xeq,  Equation  (A-43),  and 
Seq,  Equation  (A-60),  are  usually  different  for  the  two  types  of  flow.  To  remedy  this, 
the  following  definitions  for  reference  Reynolds  numbers  are  used: 


and 


where 


(A-63) 
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4  Parameters  defined  in  this  section  are  used  only  in  the  PrPr  computer  program 
and  except  for  Appendix  C  are  not  used  elsewhere  in  this  report. 
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and 
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Definitions  for  Gl,  G-p,  Jl>  xeq(  l<  and  Sgq,  l  are  given  in  preceding  sections  of 
this  appendix. 

a.  Laminar  Flow 

Using  the  definitions  given  above,  the  general  expressions  for  heat  transfer. 
Equation  (A-23),  and  skin  frietion,  Equation  (A-34),  beeome;  respectively: 
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where 
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and 


b.  Turbulent  Flow 


For  turbulent  flow,  the  form  of  Equations  (A-69)  and  (A-70)  can  be  cnanged  to  the 
Schultz-Grunow  form  given  by  Equations  (13)  and  (14).  The  final  turbulent  equations 
then  become: 
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APPENDIX  B 


SIMPLIFIED  pr-ir  METHOD 


The  complete  p  rpr  equations  presented  in  the  text  of  this  report  are  intended  for 
computer  use  only,  since  manual  computations  are  exceedingly  tedious.  However, 
simplified  equations  have  been  found  that  greatly  reduce  the  work  required  in  making 
slide  rule  calculations  with  little  loss  in  accui*acy.  First,  general  simplifications 
for  arbitrary  geometries  are  presented,  and  secondly,  equations  for  a  few  specific 
body  shapes  are  given.  Comparisons  with  results  obtained  using  the  complete  prpr 
equations  are  also  shown. 


1.  GENERAL  SIMPLIFICATIONS 
a.  Evaluation  of  prpr 


A  convenient  method  for  estimating  prpr  presented  in  Reference  2  provides 
values  within  about  three  percent  of  those  obtained  using  Equation  (A  ^12).  The 
suggested  expressions  are: 
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Note  that  Pg’Pg<  is  evaluated  at  the  stagnation  enthalpy  and  the  local  pressure, 
b.  prpr  Variatior  "  i th  x 


Unless  large  variations  in  wall  temperature  occur,  prpr  can  be  considered  to 
be  proportional  to  pressure  along  a  given  streamline.  Also,  the  effect  of  the  pres¬ 
sure  gradient  parameter  J-p  on  heat  transfer  can  usually  be  neglected.  For  example, 
the  pressure  gradient  effect  on  the  peak  turbulent  heating  rate  for  a  hemisphere  is 
only  about  3  or  4  percent.  With  these  approximations,  Equation  (8)  reduces  to: 
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and 


p  (i  =  (p  H J  (P  /P  ) 
r  r  '  r  r  m  '  e  nr 

The  subscript  m  denotes  that  the  evaluation  is  made  at  approximately  the  mean 
pressure  along  the  streamline. 

c.  Evaluation  of  E-p 


For  most  cases  the  influence  of  the  following  term  appearing  in  Equation  (32) 
can  be  neglected: 


<22c> 


exp  K 


Note  that  this  term  is  large  only  when  (x/rf)  (9rf/3x)  is  small;  consequently,  the 
influence  of  Ej  is  also  usually  small.  In  addition,  it  is  seen  in  Figure  7  that  o 
varies  only  from  .68  to  .776.  Thus,  using  an  average  or  =  .728  in  computing 


a-355  =  .8M 
r 
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in  Equation  (A -56)  will  result  in,  at  most,  an  error  of  2.4  percent.  Equation  (32) 
can  now  be  simplified  to 
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w'here 


F2,c  =  2-22^c"-294) 

Equations  (A -57)  and  (A -58)  define  2C  and  Equation  (30)  defines  Zc  Q.  Equation 
(B-5)  is  plotted  in  Figure  27  to  facilitate  computations  of  Eq-. 

d.  Evaluation  of  i.  and  Fpr 

For  real -gas  flow  in  chemical  equilibrium,  the  diffusion  influence  parameter 
can  be  obtained  either  from  Equation  (A -41)  or  Figure  20.  However,  this  term  is 
usually  quite  close  to  unity.  The  maximum  value  of  £  obtained  from  thermal  analysis 
of  the  reentry  vehicle  presented  ir  Section  vm  was  1.04,  representing  a  correction  of 
only  4%. 
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Similarly,  when  using  Equation  (20),  the  approximation  that  <tt  =  .728,  or 


.  645 
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can  lead  to,  at  most,  an  error  of  4.3  percent, 
e.  Summary  Heat  Transfer  Equation 


(B-6) 


Using  all  of  the  approximations  given  above,  Equations  (13)  and  (15)  reduce  to: 
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The  equivalent  distance  xeq  -j-  can  be  obtained  using  Equation  (B-2)  and  the 
stagnation  viscosity  from  Equation  (19). 

f.  Momentum  Thickness 


An  alternate  method  for  computing  momentum  thickness  is  provided  by: 
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If  the  effects  of  pressure  gradient  on  the  boundary  layer  profiles  are  neglected 
(Jj  =  PT  =  1.0),  the  equivalent  distances  for  heat  transfer  and  skin  friction  are 
equal.  For  this  case,  Equation  (B-9)  can  be  written: 
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2.  SPECIAL  CASES 


In  most  cases,  the  most  tedious  task  in  computing  heating  rates  is  the  computa¬ 
tion  of  the  equivalent  distance.  For  many  common  body  shapes,  Equation  (B-2)  can 
be  simplified  further.  The  following  table  summarizes  the  expressions  for  equivalent 
distance  for  some  simple  configurations. 
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Equations  for  evaluating  N  and  n  are  presented  in  Section  V  for  the  respective 
configurations. 


3.  COMPARISONS  WITH  COMPUTER  RESULTS 

Comparisons  of  Stanton  numbers  along  the  stagnation  line  of  a  sharp  ogive  at 
angle  of  attack  are  shown  in  Figure  28a.  The  simplified  PrPr  results  were  obtained 
using  Equation  (B-7)  together  with  an  equivalent  distance  computed  using  Equation 
(B-2).  These  results  are  seen  to  agree  within  2*J  with  those  obtained  using  thePrPr 
program  described  in  Volume  III.  The  free  stream  flow  conditions  represent  a 
typical  test  condition  in  the  CAL  high  energy  shock  tunnel. 


Similar  comparisons  along  the  centerline  of  a  blunt  delta  wing  are  presented  in 
Figure  28b.  The  configuration  and  flight  conditions  (V^  =  18,700  ft/sec  and  ALT  = 
187,000  ft)  correspond  to  the  orbital  reentry  analysis  described  in  Section  VI13.  The 
Stanton  numbers  obtained  from  both  methods  have  been  corrected  for  bluntness 
effects  using  Equation  (74).  The  momentum  thickness  used  in  computing  bluntness 
corrections  for  the  simplified  PrPr  values  were  obtained  using  Equation  (B-10). 
Results  from  the  two  methods  are  seen  to  differ  by  less  than  about  4%. 
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APPENDIX  C 


TURBULENT  NONSIMILAR  METHOD 


The  primary  purpose  in  developing  this  method  was  to  provide  a  means  of 
computing  complete  boundary  layer  profiles  for  compressible  turbulent  flows. 
Preliminary  results  from  this  analysis  are  reported  in  Reference  25.  The  flow 
equations  are  solved  in  partial  differential  form,  and  unlike  previous  methods, 
empirical  correlations  are  required  only  in  defining  the  Reynolds  stresses.  Solutions 
obtained  in  inis  way  are  not  subject  to  the  usual  restrictions  regarding  profile  simi¬ 
larity;  hence,  effects  of  streamwise  pressure  and  wall  temperature  gradients  are 
easily  included.  This  approach  has  been  successfully  used  in  treating  laminar  flows 
for  several  years  (References  1  and  6). 

The  turbulent  nonsimilar  equations  have  been  programmed  for  both  the  IBM  7094 
and  the  Univac  1108  digital  computers.  This  program,  including  equations  and  sam¬ 
ple  input  sheets,  is  presented  in  Volume  HI.  The  existing  program  is  restricted  to 
ideal  gases ,  but  real  gas  effects  could  be  included  by  adding  tables  for  obtaining  gas 
properties  as  a  function  of  enthalpy  and  pressure. 

The  development  of  the  turbulent  nonsimilar  method  is  presented  in  four  parts: 

1.  Derivations  of  flow  equations 

2.  Evaluation  of  turbulent  stresses 

3.  Heat  transfer  results 

4.  Convergence  and  stability 


1.  DERIVATION  OF  FLOW  EQUATIONS 

Several  formulations  of  the  fundamental  turbulent  flow  equations  can  be  found  in 
the  literature  (e.g. ,  References  3  and  7).  Normally,  the  suggested  expressions 
consist  of  the  corresponding  laminar  equation  with  an  added  term,  or  terms,  repre¬ 
senting  the  influence  of  the  turbulent  fluctuations.  The  flow  equations  can  then  be 
written  for  two  dimensional  flow  ,  in  the  following  form: 
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Tj  is  the  effective  turbulent  shear  and  is  the  effective  turbulent  heat  conduction 
term.  It  should  be  noted  that  and  q^,  do  not  truly  represent  viscous  and  conduc¬ 
tion  effects,  but  account  for  the  net  effect  of  the  turbulent  fluctuations.  Except  for 
these  turbulent  stress  terms,  Equations  (C-l)  and  (C-2)  are  identical  with  the 
corresponding  equations  for  laminar  flow. 

Minor  differences  exist  in  the  literature  concerning  both  the  definition  and 
formulation  of  the  turbulent  stress  terms.  However,  in  the  present  analysis, 
empirical  correlations  are  ustl  in  evaluating  these  terms,  and  the  exact  analytical 
formulation  is  immaterial. 

It  is  now  convenient  to  introduce  a  velocity  parameter  v  defined  by: 
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(C-4) 


T^e  streamline  divergence  parameter  r  is  introduced  into  Equation  (C-4)  to 
account  for  three-dimensional  effects.  For  laminar  flows  it  is  seen  from  the  mass 
conservation  equation  that  v  is  equal  to  the  velocity  component  normal  to  the  surface 
(v).  This  is  not  necessarily  the  case  for  turbulent  flows,  because  of  the  fluctuating 
terms  denoted  by  primes  appearing  in  the  conservation  equation  shown  below: 
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Using  Equation  (C-4),  Equations  (C-l)  and  (C-2)  can  be  written: 
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Solutions  are  to  be  obtained  by  numerically  integrating  Equations  (C-6)  and  (C-7) 
with  respect  to  x.  In  order  to  perform  this  integration,  8u/8x  and  81/9  x  must  be 
expressed  explicitly  by  the  flow  properties  at  a  given  station  in  x.  That  is,  the  right 
side  of  Equations  (C-6)  and  (C-7)  cannot  contain  derivatives  with  respect  to  x  other 
than  9P/8x  and  8r/8x,  which  are  specified  at  all  x  locations.  The  purpose  of  the 
following  derivation  is  to  obtain  an  expression  f  r  v  that  does  not  contain  x 
derivatives  other  than  8P/8x  and  8r/8x. 


Using  the  equation  of  state  (C  —3) ,  Equation  (C-4)  can  be  expanded  to: 

-  pg  -  pv  («£»  -  «£i)  -  p|U  pU  E£s)  (C-8) 

We  can  assume  8P/8y  =  0,  since  this  analysis  is  restricted  to  thin  boundary  layers. 
Also,  assuming  an  ideal  gas  (i  =  cpT) ,  Equation  (C-8)  can  be  further  expanded  to  give: 
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Making  the  indicated  cancellations  and  rearranging.  Equation  (C-12)  becomes: 


Replacing  A  and  B  in  Equation  (C-16)  using  Equation  (C-ll), 
for  V;  now  becomes: 
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the  final  expression 
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can  now  be  obtained  explicitly  from  given  u  and  H  profiles  providing  that  Tj,  q-p, 
3r/3x,  and  3P/3x  are  known.  First,  vj  is  computed  using  Equation  (C-17)  with  the 
wall  boundary  conditions  specifying  that  v0  and  (3v/3y)0  are  both  zero.  Next, 
(9v/3y)j  is  computed  using  Equation  (C-14).  This  procedure  is  repeated  until  all 
values  of  v  in  the  boundary  layer  are  obtained. 


2,  EVALUATION  OF  TURBULENT  STRESSES 
a.  Incompressible  Flow 


Published  methods  for  estimating  turbulent  shear  stresses  were  studied,  but 
were  found  to  be  unsatisfactory  for  use  in  the  present  analysis.  The  Prandtl  linear 
mixing  length  expression  provides  valid  results  only  when  constant  shear  is  assumed. 
The  correlations  used  by  von  Karman  and  van  Driest  require  that  the  boundary  layer 
be  divided  into  sublayers,  with  a  separate  correlation  for  each  sublayer.  This  pro¬ 
cedure  would  unnecessarily  complicate  the  solution  of  Equations  (C-6)  and  (C-7). 

The  general  form  of  the  correlation  for  Tj  selected  for  the  turbulent  nonsimilar 
method  is  derived  in  the  following  paragraphs. 


Assuming  incompressible  flat  plate  flow  and  that 
and  (C-6)  reduce  to: 


82u 

p  ^  <<  r^.,  Equations  (C-4) 


and 


(C  —18) 


!)H 


(C-19) 


Following  common  practice,  it  is  further  assumed  that  the  velocity  profile  and 
boundary  layer  thickness  are  given  hv 


1/7 
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(C-20) 


(C-21) 


Using  Equation  (C-20): 


du 

ay 


and 


where 
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dx 
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dx  "ld2/\dx/ 

Substituting  Equation  (C-23)  into  Equation  (C-24): 


du 

dx 


*' 1  I  d6/ax 


Also,  it  is  easily  shown  that: 


(C-22) 


(C-23) 


(C-?4) 


(C-25) 


(C  -2  (i) 


Substituting  Equations  (C-20).  (C-22),  (C-25).  and  (C-2G)  into  Equation  (C-22).  we 


rr.of  • 
fc> 


(C-27) 


Combining  terms  on  the  left  side  of  Equation  (C-27),  results  in: 


From  Equation  (C-21): 


(C-28) 


(C-29) 


Substituting  Equations  (C-21)  and  (C-29)  into  Equation  (C-28): 


Solving  for  x  in  Equation  (C-29).  we  get: 


5/4  Ke\  -S/4 

x  =  (.37)  \-y- 1  6 

Replacing  the  x  in  Equation  (C-30)  with  Equation  (C-31): 


(C-30) 


(C-31) 


(C-32) 


Equation  (C-32)  is  still  unsatisfactory  for  usage  in  the  turbulent  nonsimilar  equations, 
since  6  is  rather  arbitrarily  defined.  This  term  can  be  removed  using  the  following 
substitutions.  Differentiating  Equation  (C-22)  with  respect  to  y  gives: 


(C-33) 


Multiplying  the  right  side  of  Equation  (C-32)  by  Equation  (C-34): 


(C  —34) 


(C-35) 
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Noting  that: 


.  y 

(V  6)  =  (u/a  ) 

e 


Equation  (C-35)  can  be  expressed  by: 

dr 


(C-3G) 


where 
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We  now  h2ve  an  acceptable  equation  describing  dTj  dy.  Note  that  Equation  (C-36) 
depends  only  on  local  flow  properties  and  derivatives  in  y. 

Computations  of  skin  friction  made  using  Equation  (C-G).  with  the  expression  for 
d  Tj  /dy  given  by  Equation  (C-3G)  were  found  to  be  in  poor  agreement  with  experi¬ 
mental  data.  This  result  is  not  surprising,  considering  the  approximations  used  in 
deriving  Equation  (C-3G).  For  example,  the  contribution  of  the  laminar  shear  term 
was  neglected  in  deriving  Equation  (C-36).  but  was  included  ir.  making  skin  friction 
computations.  Also,  the  1/7  power  law  .  rofile  results  in  a  finite  value  of  du/dy  at 
the  edge  of  the  boundary  layer,  but  the  correct  boundary  condition  (du/dy)e  =  0  is 
imposed  on  the  program  equations. 

Equation  (C-3G)  is  useful,  however,  in  that  it  suggests  the  general  form  of  the 
required  shear  parameter  given  by: 


dr-] 

dy 


(C-37) 


Computer  re  Its  indicated  that  a  and  c  did  not  appreciably  influence  the  velocity 
profile.  Several  values  for  b  were  tried,  and  tne  resulting  profiles  were  compared 
with  the  empirical  results.  A  value  of  12  was  selected  on  the  basis  of  <hese  com¬ 
parisons.  Computed  velocity  profiles  for  incompressible  flow  on  a  flat  plate  using 
b  =  12  are  shown  in  Figure  29  along  with  a  1/7  power  law  profile.  The  agreement 
appears  reasonably  good,  although  some  difference  in  profile  shape  is  noted. 

The  constant  c  and  exponent  a  were  determined  on  the  basis  of  comparisons 
of  skin  friction  results  from  the  program  with  estimates  obtained  using  the  well 
established  empirical  methods.  Analysis  of  skin  friction  results  from  the  program 
indicated  that  the  variation  in  skin  friction  coefficient  with  Reynolds  number 
d(ln  Cf)/d(ln  Np  e)  is  determined  by  the  exponent  a  .  After  establishing  ’.he  correct 
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Figure  29:  COMPARISONS  OF  TURBULENT  VELOCITY  PROFILES 
FOR  INCOMPRESSIBLE  FLOW  ON  A  FLAT  PLATE 


trend  of  Cf  with  e  the  constant  c  was  selected  to  yield  the  correct  magnitude  of 
Cf.  Results  obtained  using  a  =  .  S3"  and  c  =  .054  are  shown  in  Figure  .‘10,  along  with 
the  Schultz-Grunow  (Reference  9)  <  'elation.  The  agreement  is  seen  to  be  good  for 
all  Reynolds  numbers  investigated. 


b.  Compressible  Flow 

In  order  to  extend  the  program  capability  to  include  compressible  flow  it  was 
necessary,  first,  to  determine  thermal  effects  on  the  turbulent  shear  term  9tj /3y, 
and  secondly,  to  evaluate  the  turbulent  conduction  term  qj.  The  turbulent  shear  and 
conduction  were  assumed  to  be  related  by  an  effective  turbulent  Prandtl  number 
Npr  j  in  the  same  way  as  the  corresponding  laminar  terms.  Thus, 


3i_ 

dy 


(C-3S) 


whe  re 


du  3v 


Since  the  turbulent  Prandtl  number  has  very  little  influence  on  skin  friction,  it  is 
possible  to  consider  thermal  effects  on  the  shear  stress  and  turbulent  Prandtl  number 
separately. 

Appropriate  modifications  to  the  turbulent  shear  were  determined  by  considering 
adiabatic  flow  on  a  flat  plate.  The  compressible  form  for  the  turbulent  shear  gradient 
was  assumed  to  be  of  the  form: 


The  selection  of  the  form  given  by  Equation  (C-39)  was  somewhat  intuitive,  since  no 
theoretical  basis  could  be  found  for  estimating  compressibility  effects.  The  justifi¬ 
cation  for  this  selection  must  be  based  on  comparisons  with  data. 

Program  results  were  found  to  furnish  good  agreement  with  experimental  data 
presented  in  Reference  26  with  d  =  4,  as  shown  in  Figure  31.  The  compressible  skin 
friction  coefficients  were  obtained  using  the  Schultz-Grunow  method.  The  influence 
of  Reynolds  number  on  Cf/Cf  . nc  is  known  to  be  quite  small,  and  is  neglected  in 
these  comparisons. 
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Figure  30:  COMPARISONS  OF  TURBULENT  SKIN  FRICTION  COEFFICIENTS 

FOR  INCOMPRESSIBLE  FLOW  ON  A  FLAT  PLATE 


All  experimental  data  from  Ref.  26 
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Figure  31:  MACH  NUMBER  EFFECT  ON  TURBULENT  SKIN  FRICTION 

FOR  ADIABATIC  FLOW  ON  A  FLAT  PLATE 


The  combined  effects  of  Mach  number  and  wall  cooling  are  shown  in  Figure  32. 
Qualitative  agreement  between  program  results  and  data  is  noted  regarding  trends 
with  Mach  number  and  wall  cooling,  but  predicted  values  appear  to  be  generally 
lower  than  the  data. 

Unlike  skin  friction,  heat  transfer  rates  are  strongly  influenced  by  the  turbulent 
Prandtl  number.  Since  the  turbulent  Prandtl  number  cannot  be  obtained  directly  from 
experimental  results,  the  evaluation  is  made  on  the  basis  of  the  turbulent  recovery 
factor  (r),  where: 


(C  —40) 


iaw  is  the  adiabatic  wall  enthalpy.  Ie  and  ie  are  the  total  and  static  enthalpies  eval¬ 
uated  at  the  edge  of  the  boundary  layer.  Experimental  results  (e.g. ,  Reference  27) 
have  shown  that  for  turbulent  flat  plate  flows  r  is  approximately  0.  9.  In  order  to 
establish  the  corresponding  value  of  Npr  T,  program  solutions  are  obtained  using 
NPr.T=1.0,  0.  85,  and  0.71.  The  wall  temperature  is  the  adiabatic  wall  temperature 
based  on  r  =  .  9. 


The  selection  of  Npr  T  completes  the  evaluation  of  the  Reynolds  stresses  appearing 
in  Equations  (C-6),  (C-7),  and  (017). 


Recovery  factor, 


3.  HEAT  TRANSFER  RESULTS 


Comparisons  of  heating  rates  calculated  using  the  turbulent  nonsimilar  and  the 
PrMr  programs  are  shown  below  lor  a  range  of  Mach  numbers  and  enthalpy  ratios. 


q,  Btu/ft^-sec 

M 

e 

u  ,  ft  ^sec 
e 

i  /I 
w  e 

Turbulent 

nonsimilar 

method 

PrMr 

method 

5.  06 

5946 

.502 

56.  6 

54.0 

5.03 

5923 

.  s 

13.8 

14.88 

4.77 

5840 

.204 

113.  1 

10s.  7 

1. 75 

1993 

.>02 

4.7 

4.  8 

1.74 

1989 

.501 

14.  8 

14.34 

1.70 

1926 

.207 

26.4 

23.0 

Heat  transfer  rates  from  the  two  programs  are  seen  to  agree  within  5>  except  for 
the  case  where  Me  =1.7  and  iw/Ie  =  0.207. 


Comparisons  between  turbulent  heating  estimates  from  the  nonsimilar  andprpr 
methods  for  a  hemisphere  are  presented  in  Figure  34.  The  discrepancy  seen  near 
the  stagnation  point  occurs  because  the  prpr  method  extrapolates  to  zero  at  the 
stagnation  point,  but  the  nonsimilar  method  provides  the  laminar  value  at  that  point. 
At  angles  greater  than  about  30  degrees  the  two  methods  are  in  good  agreement. 

Heat  transfer  rates  for  the  linear  wall  temperature  case  shown  in  Figure  10  were 
also  computed  using  the  nonsimilar  method.  The  nonsimilar  results  are  not  shown 
since  differences  with  the  prpr  predictions  are  so  small  that  they  would  not  be 
discernible. 


4.  ACCURACY  AND  STABILITY 

It  must  be  kept  in  mind  that  the  results  from  the  turbulent  nonsimilar  method  do 
not  represent  solutions  to  Equations  (C-6)  and  (C-7),  but  are  solutions  to  a  set  of 
algebraic  equations  approximating  the  partial  differential  equations.  First,  the  size 
of  the  grid  net  must  be  sufficiently  small  to  allow  an  accurate  description  of  the 
boundary  layer  profiles.  Secondly,  numerical  instabilities  must  be  avoided  in 
performing  the  ^tep-by-step  integration. 
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a.  Accuracy 


Solutions  obtained  near  the  leading  edge  of  a  flat  plate  are  clearly  not  valid, 
since  the  boundary  layer  is  described  by  only  a  few  points.  The  minimum  number 
of  nodes  within  the  boundary  layer  needed  to  obtain  satisfactory  solutions  can  be 
determined  by  examining  program  results  along  the  plate.  Results  shown  in 
Figure  35a  indicate  that  for  incompressible  flow  over  a  flat  plate,  skin-friction 
estimates  are  within  2  -  3fr  of  the  Karman-Schoenherr  values  when  the  velocity  at 
the  first  node  from  the  wall  is  smaller  than  75%  of  the  value  at  the  boundary  layer 
edge.  A  similar  study  was  made  for  compressible  flows  by  assuming  that  Mach 
number  and  wall  cooling  effects  on  skin  friction  are  constant  along  the  plate.  The 
computed  skin-friction  coefficients  should  then  be  proportional  to  the  corresponding 
incompressible  values  along  the  plate.  Typical  results  shown  in  Figure  35b  indicate 
that  the  accuracy  criterion  given  for  incompressible  flows  will  provide  good  results 
for  compressible  flows.  The  influence  of  the  incremental  streamwise  distance  Ax 
on  accuracy  has  been  found  to  be  negligible  provided  the  numerical  instabilities  are 
avoided. 

b.  Numerical  Instabilities 


Examination  of  computer  results  show  that  numerical  instabilities  always 
originate  near  the  edge  of  the  boundary  layer,  where  oscillations  in  velocity  and 
enthalpy  w'ere  sometimes  amplified  to  the  extent  that  the  program  results  became 
meaningless.  Although  no  rigorous  analysis  for  stability  appears  possible  at  this 
time,  an  approximate  criterion  is  presented  which  is  in  excellent  agreement  with 
observed  computer  results.  The  instabilities  are  assumed  to  arise  from  the  highest 
order  differential,  i.e.  ,  the  shear  term.  The  iaminar  shear  contribution  is  assumed 
to  be  negligible  compared  to  the  turbulent  value.  The  momentum  equation  at  the 
boundary  layer  edge  then  becomes 
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(C-45) 


Initially,  the  velocity  at  one  node  is  given  by  u  =  uc  t  where  <  is  the  initial 
error  in  computed  velocity  at  y  =  yc  . 
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Figure  35:  EFFECT  OF  Ay  ON  PREDICTED  SKIN-FRICTION  COEFFICIENTS 
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N  is  now  defined  as  the  number  of  nodes  from  the  wall  to  y  .  Hence, 


(C-46) 


Neutral  stability  occurs  when  the  difference  between  u  at  y  and  the  two  adjacent 
nodes  remains  unchanged  at  the  adjacent  downstream  station,  then  the  velocity 
gradient  at  the  adjacent  nodes  is  given  by: 
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(C  —47) 


then 


The  velocity  profiles  at  the  two  x  locations  are  sketched  below: 


(C  -48) 


It  is  seen  that  the  difference  between  uv 

y  c 

locations  when 


and  u 


yt  +  Ay 


is  constant  at  the  two  x 


(C-49) 


where  Axns  is  the  Ax  required  for  neutral  stability.  Substituting  Equation  (C-49) 
into  Equation  (C-46)  gives: 


(C-50) 


or 


_  12.26 

Xns  ...833 
N 


(C-51) 


Thus,  it  is  concluded  that  if  Ax  exceeds  the  value  given  by  Equation  (C-50)  small 
oscillations  in  velocity  will  be  amplified.  Since  small  oscillations  are  unavoidable 
because  of  round-off  errors,  it  is  concluded  that  stability  can  be  maintained  only  if 
Ax  <  Axns.  This  conclusion  is  supported  by  computer  results  for  two  cases  shown  in 
Figure  36.  The  points  presented  as  unstable  are  from  the  first  x  location  at  which 
oscillations  are  observed. 


Equation  (C-51)  defines  the  maximum  incremental  distance  Ax  for  neutral 
stability.  In  specifying  the  grid  size  to  be  used  for  normal  operation  of  the  turbulent 
nonsimilar  program  it  is  recommended  that 


Ax 
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(C-52) 


where  Nmax  is  the  number  of  points  to  be  computed  at  the  most  downstream  x 
location  to  be  considered. 
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Figure  36:  NUMERICAL  STABILITY  CRITERIA  FOR  NONSIMILAR 

BOUNDARY  LAYER  PROGRAM 
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